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Abstract  of  Dissertation  Presented  to  the  Graduate  Council 
in  Partial  Fulfillment  of  the  Requirements  for  the 
Degree  of  Doctor  of  Philosophy 


A LARGE  DEFLECTION  THEORY  OF  ORTHOTROPIC 
SANDWICH  SHALLOW  SHELLS 

By 

Nai-Chi  Wu 
April,  1965 

Chairman:  Dr.  I.K.  Ebcioglu 

Major  Department:  Engineering  Science  and  Mechanics 

A theoretical  derivation  of  the  equilibrium  equa- 
tions for  a large  deflection  of  sandwich  shallow  shells 
under  mechanical  as  well  as  thermal  loading  has  been 
carried  out.  The  approach  is  in  general  coordinates 
so  that  it  is  inclusive  enough  to  cover  any  specific 
reference  frame  for  the  actual  engineering  problem.  By 
using  the  principle  of  minimum  potential  energy  five 
coupled  nonlinear  differential  equations  in  tensorial 
form  are  obtained. 

In  order  to  demonstrate  the  applicability  of  the 
above  theory,  a typical  practical  problem  with  axisymmetry 
is  treated  numerically.  By  using  a perturbation  method 
these  nonlinear  differential  equations  are  transformed 
into  several  sets  of  linear  differential  equations.  With 
the  aid  of  power-series  expansion  technique  the  equations 
are  solved  successively.  Finally,  the  problem  is  reduced 


x 


to  a simpler  case  the  numerical  solution  of  which  is 
available  for  comparison.  The  results  turn  out  to  be 
in  good  agreement  with  the  known  solutions. 
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INTRODTJCTION 


It  has  been  known  for  a long  time  that  sandwich 
construction  is  an  efficient  means  of  obtaining  a light- 
weight structure  with  high  strength.  Modern  developments 
in  the  aerospace  industry  have  created  an  urgent  need  for 
such  construction.  Economical  manufacturing  methods  are  now 
being  developed.  As  a result,  more  research  work  concerning 
these  structures  has  become  increasingly  desirable. 

por  the  sandwich  shell  the  exact  solution  of  the 
problem  is  very  difficult  to  obtain  because  in  this  case 
the  core,  as  well  as  the  facings,  has  to  be  treated  as  a 
three-dimensional  body.  Numerous  approximate  theories  in 
line  with  ordinary  shell  theory  have  been  proposed  in  the 
past.  Among  these  the  first  significant  contribution  to 
an  understanding  of  the  behavior  of  sandwich  shells  was 
made  by  E.  Reissner  (1),  who  assumed  that  the  facings  do 
not  sustain  any  bending  moment,  i.e.,  the  in-plane  stresses 
are  evenly  distributed  across  the  thickness  of  the  facing, 
and  the  core  is  soft.  The  transverse  compression  and  shear 
deformation  of  the  core  are  considered.  In  the  same  paper, 
these  effects,  which  differentiate  the  sandwich  shell  theory 
from  the  ordinary  shell  theory  based  on  the  Kirchhof f-Love 
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hypothesis , were  evaluated  for  several  types  of  sandwich 
structures.  It  was  found  that  for  certain  geometrical 
configurations  and  material  properties  the  effect  of  trans- 
verse normal  stress  deformation  is  negligibly  small.  Within 
this  approach,  Reissner  also  initiated  a finite-deflection 
theory  of  sandwich  plates  (2,  3) , based  on  the  assumption 
that  the  core  is  in  the  small-deflection  region  and  the 
facings  act  as  membranes  experiencing  a finite  deflection. 
The  resultant  equations  permit  the  analysis  of  the  effects 
of  transverse  shear  deformation  and  the  transverse  normal 
stress  deformation  in  the  core  on  the  over-all  behavior  of 
the  sandwich  plate.  It  was  shown  through  an  order-of- 
magnitude  analysis  that  the  effect  of  transverse  normal 
stress  in  the  core  is  negligibly  small  compared  with  the 
effect  due  to  transverse  shear  stresses.  It  was  further 
shown  that  the  range  of  deflections  for  which  the  linear 
theory  (small-deflection  theory)  is  adequate  decreases 
according  to  an  explicit  formula  as  the  cere  is  made  softer 
relative  to  the  facings.  Since  then,  numerous  papers  con- 
cerning the  structural  properties  of  sandwich  plates  have 
been  published.  Some  significant  contributions  (4,  5,  6, 

7,  8,  9)  have  been  summarized  in  (1C).  Wang  (U)  proposed 
a simplified  method  to  derive  the  stress-displacement 
relation  for  the  finite-deflection  theory  by  using  the 
principle  of  complementary  energy.  The  result  was  applied 
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to  set  up  the  governing  equations  for  cylindrical  sandwich 
shells  and  also  to  study  the  stability  of  sandwich  cylin- 
ders {12,  1J3) . A rather  complete  survey  of  modern  develop- 
ments in  the  analysis  of  sandwich  structures  has  been  con- 
ducted by  Habip  (14).  In  the  field  of  large-deflection 
theory  of  sandwich  shallow  shells,  Karavanov  (15)  derived, 

by  following  Reissner°s  approach  (2,  3 ) , the  nonlinear 

deflection  equations  taking  account  of  the  curvature  effects. 
Grigoliuk  has  derived  a set  of  governing  differential  equa- 
tions for  sandwich  shells  with  a soft  core  (_16) , by  con- 
sidering that  the  facings  transmit  both  the  moments  and 
tangential  stresses  in  the  finite-deflection  region  and  the 
core  transmits  only  the  transverse  shear  in  the  small- 
deflection  region.  Later  he  also  considered  a case  of 
finite  deflection  of  sandwich  shells  with  a stiff  core,  in 
which  the  membrane  stresses  were  also  taken  into  account 
(l_7 ) . The  same  approach  has  been  worked  out  by  Wasicina  ( 18 ) 
for  the  case  of  unsymmetrical  construction.  In  1961,  Mush- 
tari  (19)  published  a paper  to  discuss  a general  theory  of 
shallow  sandwich  shells.  In  this  paper,  an  unsymmetrical 
sandwich  shell  construction  was  proposed.  The  thermal 
effect  was  also  taken  into  account.  However,  the  temperature 
was  limited  to  uniform  distribution  over  the  facings  and 
linear  distribution  through  the  thickness  of  the  core.  In 
the  same  year,  Fulton  (20)  presented  a paper  on  a nonlinear 
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theory  of  shallow  unsymmetrical  sandwich  shells  of  double 
curvature.  Grigoliuk  and  Chulkov  have  presented  a small” 
deflection  theory  of  sandwich  shells  with  stiff  core  (21 ) . 
They  also  considered  the  compressibility  in  the  transverse 
direction  of  the  core.  The  core  was  taken  as  a three= 
dimensional  body  but  with  the  assumption  that  the  displace- 
ments could  be  approximately  expressed  as  a linear  function 
of  the  transverse  coordinate.  Both  statical  and  dynamical 
cases  were  discussed.  However,  only  the  governing  equations 
were  set  up  and  no  solution  to  the  equations  was  made  avail- 
able. Grigoliuk  and  Chulkov  have  published  a paper  on  a 
general  large-deflection  theory  of  elastic  sandwich  shallow 
shells  (22).  In  contrast  to  the  paper  given  before,  this 
paper  dealt  with  an  elastic  shallow  shell  with  different 
thicknesses  of  the  facings.  The  facings  were  considered 
homogeneous  and  isotropic  but  of  different  materials.  The 
core  was  assumed  incompressible  in  the  transverse  direction, 
and  the  displacements  in  the  tangential  directions  were 
assumed  to  be  a linear  function  of  the  transverse  coordi- 
nate. The  reference  surface  was  defined  in  terms  of 
rectangular  Cartesian  coordiantes.  Thus  far,  the  Russian 
works  mentioned  above  do  not  include  any  solution  to  the 
governing  equations. 

This  dissertation  deals  with  a large-deflection 
theory  of  sandwich  shallow  shells  under  mechanical  as  well 
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as  thermal  loadings.  The  analysis  is  carried  out  in  general 
coordinates.  All  the  parameters  are  in  terms  of  functions 
of  the  latter.  Naturally,  tensor  notation  is  employed  in  the 
derivation.  Through  differential  geometry  the  strains  can 
be  expressed  in  terms  of  the  displacements,  including  the 
nonlinear  terms,  of  the  middle  surface  of  a shell  within 
an  approximate  theory  of  shallow  shells.  This  relationship 
is  applied  to  calculate  the  strain  energies  for  the  upper 
and  lower  facings  as  well  as  for  the  core  which  is  con- 
sidered as  stiff  and  orthotropic.  The  total  strain  energy 
is  obtained  by  summing  those  strain  energies  for  the  in- 
dividual components  of  the  sandwich  structure.  The  varia- 
tion principle  of  minimum  potential  energy  is  employed  to 
obtain  a set  of  nonlinear  equations.  Certain  assumptions 
used  in  this  derivation  are:  (a)  the  facings  and  core  are 

thin,  and  (b)  the  Kirchhof f-Love  hypothesis  is  valid  for 
the  facings  while  the  tangential  displacements  for  the  core 
are  linear  functions  of  the  thickness  coordinate.  The 
transverse  deflection  is  assumed  constant  through  the 
thickness  of  the  whole  structure. 

Two  methods  are  proposed  in  the  derivation.  In 
the  first  method,  labeled  method  (A),  the  tangential  dis- 
placements of  the  middle  surfaces,  u^  , u£  , of  the  upper 
and  lower  facings,  respectively,  and  the  transverse  deflec- 
tion w are  taken  as  the  independent  parameters,  as  most 
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previous  works  have  done.  The  second  method,  labeled  method 
(B),  takes  , the  tangential  displacements  of  the  middle 
surface  of  the  core,  , the  rotations  of  the  cross-section 
normal  to  the  middle  surface  of  the  core,  and  w,  the  trans- 
verse  deflection  of  the  structure,  as  the  independent  para- 
meters. Both  methods  have  specific  merits  in  practical 
applications.  Since  there  is  no  difference  in  principle, 
only  method  (b)  is  applied  in  a practical  application  in 
which  all  the  tensorial  equations  are  reduced  to  engineering 
form. 

An  axisymmetrical  sandwich  shallow  shell  of  para- 
boloid  geometry  is  introduced  as  an  illustrative  example  to 
demonstrate  the  applicability  of  the  proposed  theory.  A 
perturbation  scheme  using  the  dimensionless  central  deflec- 
tion as  a parameter  is  chosen  to  attack  the  nonlinear 
equilibrium  equations  obtained  from  the  previous  theoretical 
derivation.  With  this  scheme  the  nonlinear  differential 
equations  are  reduced  to  several  sets  of  linear  differential 
equations.  In  solving  these  linear  differential  equations 
certain  series  solutions  with  undefined  coefficients  are 
successively  substituted  into  the  differential  equations. 
Together  with  the  boundary  conditions,  a set  of  linear 
algebraic  equations  of  the  undefined  coefficients  are  ob- 
tained. The  matrice  inversion  is  applied  to  solve  for  the 
undefined  coefficients. 
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A digital  computer  program  written  in  Fortran  II 
provides  the  numerical  results.  Furthermore,  the  case  is 
reduced  to  an  ordinary  shallow  spherical  shell  to  make 
some  numerical  comparisons  with  known  solutions  (23,  24,  25). 
xhis  leads  to  good  agreement. 


CHAPTER  I 

THEORETICAL  DERIVATION 
1 . Mathematical  Preliminaries 

The  geometrical  relations  in  plane  coordinates 
for  a shallow  surface  have  been  discussed  by  Green  and 
Zerna  (26).  A net  of  coordinate  curves  on  the  middle 
surface  M is  considered  to  be  projected  on  a plane  K. 

The  direction  of  projection  is  given  by  a unit  vectos  e3 
which  is  perpendicular  to  the  plane  7T . Hence,  there  now 
is  a one-to-one  corresponding  net  of  curvilinear  coordinates 
between  the  7T-plane  and  the  surface  M.  The  position  vector 
P can  be  expressed  as  shown  in  Figure  1, 

R = 7 * zes.  (1* 

In  this  case  the  surface  M can  be  represented  by  a net  of 
curvilinear  coordinates  on  the  7T“plane  and  a scalar  quan- 
tity Z which  is  a function  of  the  curvilinear  coordinates 
on  the  7f=plane. 

The  tangent  vectors  to  the  curves  on  the  surface  M 
can  be  expressed  by 

a.,  = - € ' -C-  ' (1* 
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where  r 7,  and  the  comma  means  partial  derivative 
with  respect  to  the  indicated  variable. 

The  displacement  vector  can  also  be  expressed  by 


U - <X  + w a ; - e*  v-  A/  <ff  , c<=/.  2 , (1-3) 

where  a*  , i = 1,2,3,  and  e*  , i - 1,2,3,  are  the  right 
hand  triads  on  the  surface  M and  the  7T-plane,  respectively, 
and  , w and  u^  , w are  the  displacement  components  with 

respect  to  the  base  vectors  "a*  and  e^,  respectively.  For  a 
shallow  shell  it  is  understood  that  the  values  of  Z and  its 
derivatives  up  to  the  second  order  are  small  in  comparison 
to  the  other  dimensions  of  the  shell.  Then,  as  an  approxi- 
mation, we  neglect  terms  containing  the  squares  and  products 
of  Z as  well  as  its  derivatives  up  to  the  second  order  with 
the  other  quantities.  Thus, 

£2  = ( Q ) ( ^3  * ^ ) ~ <•  ( 1 -^ ) 

Similarly, 

cf?  ± e*?  (1-5) 


and 


a.  = 


(1=6) 


1C 


Multiplying  both  sides  of  Equation  (1-2)  by  e°^  and  applying 
Equation  (1-4),  we  can  obtain 


y- 


(1-7) 


From  vector  algebra  we  have 

aa  x as  = ± a3  , (1-8 

X • (1'8 

Substituting  Equation  (1-2)  into  Equation  (1-8) , and  assigning 
o<  = 1,  /3  = 2,  we  have 


— ( &,  + <^3  ) *(  * Z,2  € ) . 


(1-9) 


Therefore, 

^ 3 ~~  ^3  ~ ' <3  (~>  » 

or 

— ^3  o/  ^ 

Using  Equation  (1-7),  Equation  (l-ll)  becomes 

~ + Z,  ^ 5 ) • 


(1-10) 


(1-11) 


(l-12a) 
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Neglecting  terms  involving  second-order  products  of  Z-terms, 
we  obtain 

e3  = a3  + ZJO,  a*  . (1-I2b) 

In  order  to  calculate  the  strain  components  in 
terms  of  the  displacement  vector  we  define  the  displacement 
vector  (see  Figure  2)  as 

)/  - A - T / . (1-13) 

The  Lagrangian  strain  tensor  is  defined  as  (27) 

= T( <?.j  - #.-j)  . (1-14) 

Where  Gjj  and  g^j  are  metric  tensors  of  the  initial  and 
deformed  body,  respectively.  Performing  the  partial 
differentiation  in  Equation  (1-13)  we  obtain 


Also,  by  definition 


- ’ % + /v  * Vj  * ^ ' ^7  • 


(1-15) 


(1-16) 


Substituting  Equation  (1-16)  into  Equation  (1-14),  we  ob- 
tain the  Lagrangian  strain  tensor 
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s‘j  = * ^ Z * KyZj), 

and  on  a surface  the  strain  tensor  can  be  expressed  by 


(l=17a) 


^ ~ 2 ( • l<j3  + ap-V'o,  + V*'  )• 


(l-17b) 


Referring  to  the  surface  M,  we  can  also  write  the  displace- 
ment vector  in  the  form 


V = U*  a + n as  , 

and  referring  to  the  7T-plane,  as 


( l-18a ) 


• o/ 


1/  = u*  e + u . 


(l-18b) 


Hence 

Un  / d dj  = e*  + A/  e3  . ( l-19a) 

Using  Equation  (1-7)  and  Equation  (1-11),  we  find 


a*  + d a3  = (u^  + e^)d3.{l’l9b) 

Hence 

U = U.  i-  h/  Z 


(l-20a) 
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IaJ  = W ~ U*  . (l-20b) 

f 

Also 

~ ua/^  e * / M,*,  <?3  > (1-21) 

where  (/)  means  covariant  derivatives  with  respect  to  the 
indicated  variable.  Substituting  Equations  (1=2)  and  (1-21) 
into  Equation  (1— 17b) , we  obtain  the  strain  tensor 


_ / 


T3 


U. 


Since  the  displacements 


+ Uf/o<  + 2 'ft 

in  the  middle  surface 


are  still  small 


even  though  we  consider  a large-deflection  case,  by  neglect- 
ing the  higher  order  terms  u ^ u///3  e in  comparison  with 
the  other  terms,  we  can  express  the  strain  tensor  in  terms 
of  the  displacements 


orresponding  strain-displacement  relations  have  been  given 
to  rectangular  Cartesian  coordinates  by  Marguerre  (28),  and 
in  cylindrical  coordinates  by  Nash  (23).  Here  an  alternate 
derivation  has  been  presented. 

2.  Derivation  of  Equations 
A.  Method  A 

xhe  displacement  of  the  middle  surface  of  the  upper 
facing,  u8^  , and  the  displacement  of  the  middle  surface  of 
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the  lower  facing,  u'^  , and  the  transverse  deflection  w are 
taken  as  the  reference  variables. 

a.  Strain-Displacement  Relations  and  Strain  Energy 
of  a Sandwich  Structure  (Figure  3) 

The  facings  are  considered  thin  and  the  Kirchhoff- 
Love  hypothesis  is  valid.  The  displacements  for  the  facings 
can  be  expressed  by 


O 


(1-23) 


where  are  the  displacements  at  the  middle  surface  and 
Xj  are  the  coordinates  of  the  middle  surfaces  of  the  facings. 
For  the  upper  facing,  with  u^  = u^  , x§  = (t  + t),  and 
^ t ^ x3^t  + the  displacements  can  be  written  as 


u'  ~[ *3  - z (t  + t)]lV, 


o 


(1-24) 


The  strains  can  be  obtained  by  substituting  Equation  (1-24) 
into  Equation  (l-22b) 


- fe  - £(t  / t)]k/A 


r 


(1-25) 


where 


^ = i( u-l/p  + i-  z.„ Up  t / k,  fty) . (1'26) 

For  the  lower  facing,  with  u°,  = u’^  , ^ (t  + t),  and 
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^ 2 ^x3  2’  disPiacemen'ts  can  be  written  as 

-[ x3 +£(*  + *>] w <*  • u-27) 

The  strains  can  be  obtained  by  substituting  Equation  (1-27) 
into  Equation  (l-22b) 

S<*p  = -A  +±(t+i , <1-28) 

where 


si  = ^ 


7* 


a. 


*//* 


/ 


( 1-29) 


The  variation  of  strain  energy  for  the  upper  facing  is 

/ U ds  =/ -Uj  - 

Where  cr/7J  is  the  stress  tensor  for  the  upper  facing, 
he  variation  of  strain  energy  for  the  lower  facing  is 


su'd  -(z+H&s*Jsk  «i-m) 


y3 


r 


rj 


where  cr"7*  is  the  stress  tensor  for  the  lower  facing. 

f-or  the  core,  — ^ < x3  ^ and  since  we  consider 
a stiff  orthotropic  core,  we  have  to  consider  the  strain 
energy  due  to  the  deformation  of  the  middle  surface  and  the 
transverse  shear  deformation.  Under  this  circumstance,  the 
Kirchhof f-Love  hypothesis  no  longer  holds  for  the  core. 
However,  it  is  valid  to  assume  that  planes  normal  to  the 
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middle  surface  of  the  core  remain  plane  after  deformation  but 
no  longer  necessarily  normal  to  the  middle  surface  because 
of  transverse  shear  deformation  in  the  core.  Under  this 
assumption,  the  displacements  in  the  core  are  linear  func- 
tions of  the  distance  measured  from  the  middle  surface  of  the 
core.  We  let  u ^ = A + Bx^,  where  A and  B are  constants  to 
be  determined  by  applying  the  geometrical  continuity  at  the 
interfaces.  At  the  interface  between  the  upper  facing  and 
the  core,  i.e.,  at  *3  = £•  , the  displacements  are  given  by 
Equation  (1-24).  Thus 

A f 8 Sr  = . (l-32a) 

At  the  interface  between  the  core  and  the  lower  facing,  i.e., 
*3  = ” |,  from  Equation  (1-27)  we  have 

A - B i = U"  ~ ~ W/oi  . (l-32b) 

Solving  Equations  (l-32a)  and  ( 1- 32b ) for  A and  B,  we  obtain 

A - ( UL<  + u"  ) , (l-33a) 

3=  T(U^  ~ U'~  + (l-33b) 

Hence,  the  displacements  in  the  core  can  be  expressed  as 

u*  = z(u~  + UD  + y~( - u*  * . 


(1-34) 


17 


The  strains  can  be  obtained  by  substituting  Equation  (1-34) 
into  Equation  (l~22b)  so  that 


£<*/3  - z[ -2  ( u-/p  ? u-y/, 


+ i(*. 


'fc,  r u ,*/<*)  y f'(  ur^  y tJ/^) 


+ z/a,  w/p  + £/*  A4,  y 


'/*  yr 


(1-35) 


The  transverse  shear  strains  can  be  obtained  from  Equation 
(1-17)  as  follows 

Sas  = i(  + %,■$,).  (1'36) 

Since  V = u^  "a^  + wa^,  we  have  V,  3 = u ^3  a”*  + w,3  T3; 

the  last  term  equals  zero  since  w is  assumed  to  be  indepen- 
dent of  X3.  So  that  we  have 

Kj  = «V,  aw  . (1-37) 

Substituting  Equations  (1-21)  and  (1-37)  into  Equation 
(1-36),  we  obtain  the  transverse  shear  strains  in  terms  of 
the  displacements 

cFy)  . (1-38) 

Dy  the  same  argument  utilized  before,  the  higher  order  terms 
u ^ u y/ 3 a^y  are  neglected.  Then, 
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£ & 3 ■<?  ( ^°</3  ^ at  ) ■ 


(1-39) 


When  Equation  (1-34)  is  used,  Equation  (1-39)  becomes 

<£.,  +(t+  i)^] . 

.-'.fv-  . 

The  variation  of  strain  energy  for  the  core  is 


r 


S u ^ Jl/ 

=/ +«T.)  + 
y + ^P//3f  *£< 

1 1 u~  - ->- ( t + iJ  W/a  ]3l/  , 


-43 

cr 


(1-40) 


(1-41) 


where  P^/3  are  the  normal  stresses  and  tra'3  are  the 
transverse  shearing  stresses  in  the  core. 

b . Derivation  of  the  Equilibrium  Equations 

The  theorem  of  minimum  potential  energy  (29)  is 
applied  to  derive  a system  of  equilibrium  equations. 

The  total  potential  energy  for  the  sandwich  shallow 

shell  is 

p£  = u'+u"+U - 


(1-42) 


(9fr-T) 


?p  i f - 


sp  $ 2>"  +("'?* 


pZ 


>/^ / ^4  rr  {Mf  + 


m?  ^[r/? 4 4 >/  t-'^JJ-  =-n? 

sauiooaq 

(^-1)  uofqenbn  ‘paxxdcle  sx  uiajoaqq  aoua6.iaAXP  aqq  uaqv. 


(qQW) 


• r'r(jh-r')+~ 


r 

1 


1 + 


(eQfr-I) 


7f=,- 


■-*  j=j~-N 
f -f-/ 


a jaqM 


(trt'-T) 


'rpfaf&v-  C^s 4 y" 4 J = 

'ir[‘//%’f(jrj  ■ r'--  4 c^s’/p*  *'/~r , /?  ? 


(oc-l)  uoTqcnba  uioaj 

(Efr-T)  ‘ o = f-  OS  f + M $ = ("3  d )f 


saqeqs  A6jaua  xeT^ue^oc^  uunujiuxui  jo  wajoaqq  aqj. 


61 
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Similarly  from  Equation  (1-31) 

su"= 

+ Artif  Jhtj ja  * fs  [ /V" 

/ /V"“%£  i-z/Ji-  /*rtj  ]S»J yj 
-£  Ari’j*t£i/,ds. 


where 


_ t 

! * 
<*/i  ! 


cr' 


^ g/Sj  , 


_ / 


r* 


M"0!3  = i 


/_  t 


J 


ft 


(1-47) 


( ] -48a ) 


(l-48b) 


and  from  Equation  (1-41),  we  have 

iU  = -lf/v^  ±(fy.  * /<;  //>^ 

+ T a ^ ( /<  - S“'L)- 

- f - /<;  ' T1^  t”}**  11 

+ £ f /v^fifSK  +S^)t(^+>&)Sw]  + F 
(SFV  - Fa")  - /ax  Sk/  f ^ ^ 7~ V? 
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where 


t 

2 

£ 


(l-50a) 


(l~50b) 


, — i ^ -3  j — o<  J 

Q = I . ^ c/*9  . 

t 3 


(l-50c) 


Substituting  Equations  (1-46),  (1-47)  and  (1-49)  into 
Equation  (1-43),  we  obtain 


-/ {(^  ' - t<r")Sa. 

* ( A/"  Cja  + ^ W ^ ~ T ^ 'jj  + T 9 £ < 

+ If  ( Ar°f+  ""f*  * r'X  *%,  * uif  + ( /T^ 

t 


/ 

C< 


- 


t t t £V  3 

? # ify  / J / 


v> 


+f  { f(A/'^+  t ( #*-9*+ jr? 


~ ? +f  ( AT*/3  + /V'r//J+  /V  7*J( + ^)+(Af'Xz, 

tsr'°^  - t 

/ §~  j jjz  Js  - o . 


(1-51) 
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Since  Ju9  , /u'^  and  /w  are  independent  variations 

the  coefficients  of  these  terms  in  the  surface  integral 


should  vanish  in  order  to  satisfy  Equation  (1=51). 
we  have  the  five  simultaneous  equilibrium  equations 

Hence , 

{/v^  + */v”'/2)sfi  + -o°'3)=0. 

(l-52a) 

(a/"^3  + ^ ~ ~ Q°'S')  = ° ' 

(l-52b) 

+ ?<?*+ 7-0  , 

where 

(l-52c ) 

N*<a  = tf'V  + sr'l’ + , 

( 1 -53a ) 

M*?  = t m"°*‘  / /if~P  , 

(l-53b) 

? tr  £ 

* = 7 ■ 

(l-53c) 

and  the  boundary  conditions 

(/V'^  + i Jp  =o 


(l-54a) 
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(AT?  f mV) Sul  £ = o , 

/ /T* r 4,  A ^ - r/V^  * t(fs)]^  Sn/  = o , 

M'^  -f  A1"°'/'1  -/■  f At’1'3)  J "*£,  A = o. 

So  far,  we  have  obtained  a set  of  equilibrium  equa- 
tions and  boundary  conditions  in  terms  of  the  stress  and 
moment  resultants.  Now  we  want  to  express  these  equations 
and  boundary  conditions  in  terms  of  the  displacements.  Let 
us  use  a generalized  form  of  the  Duhamel-Neumann  equations 
for  an  orthotropic  material  under  the  generalized  plane- 
stress  assumption 

o-V  = /T'/’ T . (l-55a) 

0-“*  = , (l-55b) 

where  kV'-1  and  k‘'J/J  are  the  elastic  coefficients 
which  depend  on  the  metric  tensor  as  well  as  the  properties 
of  the  undeformed  body.  T is  the  temperature  function. 

D°^  are  thermal  coefficients.  Substituting  Equations 
(1-25),  (1-26)  and  (1-55)  into  Equations  (l-45a)  and  (l-45b), 


(l-54c) 

(l-54d) 
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and  using  the  symetric  properties  of  the  coefficients,  we 
have 


u>/j  f i'f  % +£ty,>V/s) 

-(*3-  u-56a) 

= /tw'3rsft( u'p  * £ ^ r^, 


where 


J+t 


T'  ^ 


2 


M 


■ */3 


+ tr 

j at/SYjf  f 
/ ^ ' /(UVj  * z4  W/j 


(l-56b) 


(*3  ~ t f*3  - 'p-)^S3 


(l-56c) 


- _ , Ajt.  of/a 

~ * /<?  WyS  / Afr  ' J 

~ S'  i ' t 

M'r  =J£  if'/9 T' ( sT3  - JsT3  . (l-56d) 

Similarly,  by  Equations  (l-48a)  and  (l-48b),  we  have 

/V"  ''  = ^ ( U>%  ^4y^  + * + T"*/3 \ (l-57a) 
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where 


7" 


_ i 


P*/3  r'J*,  , 


A/' 


r=  - 4%,  / #■?  , 


-I 

*2 


A/"  ' =/  £>°'/3t*(. ^ y-  ArA 


Also,  applying  Equations  (l-55a),  (l-55b)  and  (1-40) 
Equations  (l-50a),  (l-50b)  and  (l-50c),  we  obtain 


t 

— cS/2  [ / 

/V  7 =/ 

Lt  t 


+ t % ) -f  ± *%,  w h tv/jj* 

- A~'*^ tj  jr(  u'p.  / U'^)r  Jr tv/rH,/s  +%,!%]+ 


where 


£ 

r 

T ' - / _ J?  7 7*  c/X, , 
J - t ’ 


*/*  - t V ,/ 


a * t^rs)  t ' 


(l-57b) 


(l-57c) 


(l-57d) 

to 


(l-58a) 

rr. 


(l-58b) 


( l-58c ) 
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t 


2 


D*'3  T 


and 


(l-58d) 


Q™3  - J * l u>' a>  +( t + t)  ¥$>] . (l-58e) 


From  Equation  (1-53),  it  can  be  written  that 


/v"'’  = A*/* rJ f t ( u'vj  + u'Vj  +2  2^*4  + /*}.*%)] 
*/)  7 t f j (u-fy  * u"^)  -t  £ 

^ z4*/j]  1-  T'^ / 7 71 ' 

/V^=  -/> / A *“£(«*  - 

-+  ) / AtCf  + + At 


(1-59) 


(1-60) 


By  substituting  Equations  (1-56),  (1=57),  (1-58),  (1-59)  and 
(1-60)  into  Equations  (1-52)  and  (1-54),  we  obtain  five 
simultaneous  differential  equations  and  boundary  conditions 
in  terms  of  the  displacements  u ^ , u'^  and  w. 


B.  Method  B 


Since  the  equilibrium  equations  obtained  from  method 
(A)  have  a rather  complicated  form  in  terms  of  the  reference 
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displacements  , u"^  and  w,  it  is  desirable  to  convert 

these  equations  into  a simpler  form  in  order  to  solve  the 
equilibrium  equations.  Different  reference  variables  are 
chosen  here  to  derive  the  equilibrium  equations. 

We  take  u^  , the  displacement  of  the  middle  surface 
of  the  core,  ©<*  , the  angular  displacement  of  the  cross 

section  originally  normal  to  the  middle  surface  of  the  core, 
and  w,  the  transverse  deflection  of  the  sandwich  shell,  as 
the  reference  variables. 

a.  Strain-Displacement  Relations  and  Strain  Energy 
of  a Sandwich  Structure 

^he  same  assumption  made  in  method  (A)  is  used 
for  the  core;  i.e.,  planes  normal  to  the  middle  surface 
of  the  core  still  remain  plane  after  deformation.  Consequently 
we  have  the  displacements  for  the  core  expressed  by 


at 


= 

+ *3 

>> 

(1-61) 

^ = 

t 

2 ; 

= 

(l-62a) 

^ = - 

t 

-?  ; 

^ = 

t 

2 

4,  . 

(l-62b) 

Substituting  Equation (1-61 ) into  Equation  (l-22b). 


we  can  express  the  strains  in  terms  of  the  displacements 
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^ (l‘63) 
y equation  (1-39),  we  obtain  the  transverse  shear  strains 

= i ( &<*  + . (1-64) 

he  variation  of  strain  energy  for  the  core  can  be  obtained 
as 

</  U = f V/'  / 2^ JV 

=i + iycwf  (1-65> 

* t - t v&jJj a , 

where  1 , M ^ and  Q0^3  are  defined  by  Equations 

(l-50a),  (l-50b)  and  (l-50c),  respectively. 

For  the  upper  facing,  with  | < x3  ^ | + t , and 
applying  the  Kirchhof f-Love  hypothesis,  we  have  the  dis- 
placements expressed  as 

U^-U^  + / 0^  - (x3  - /;  ^ . (1-66) 

From,  Equation  (l-22b)  the  strains  in  the  upper 
facing  can  be  found  in  terms  of  the  displacements 
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+ +tyHu  *&*$,]. 


(1-67) 


For  the  lower  facing,  with  - -~t^'xo<C-  1 

2 J v 2 

and  applying  the  ^irchhof f -Love  hypothesis,  we  have  the 
displacements  expressed  as 


~ ^cK  2 ~ ( 1"  £ ) Ifytx 


(1-68) 


For  the  strain-displacement  relations  of  the  lower 
facing,  using  Equations  (1-22)  and  (1-68),  we  obtain 


<4^  - z [ U°f  a 0o.//3' 2(Z,+  y)^  + Up/* 

2 @/3/<x  + % + */&  Wfr  + ^Ux  Wf?  J . 


(1-69) 


is 


The  variation  of  strain  energy  for  the  upper  facing 


t>U'=jr 


^ ^ L * 2 &<*/a  ('tj 


7* 


°r 


7 j /V' / U f4&  % / * 


(1-70) 


+2  0^ 


where 
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t 


(l-71a) 


= /_  . 

'he  variation  of  strain  energy  for  the  lower  facing  is 

S U"  =fA  f tr'P  S( 

* i ^ Jc/a  , 


ftt 


(l-71b) 


(1-72) 


where 


/V" 


-t 

*2 

'-I-* 


•VS*  > 


(l-73a) 


. t 


/V' 


t 


<7^(4 1 i)  ^ 


(l-73b) 


b . Derivation  of  the  Equilibrium  Equations 

The  total  variation  of  the  strain  energy  stored  in 
the  sandwich  shell  is 
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<n/r  = ji/'+s  l/''+jd 

~JA  / ' ^ 7 74  74  2 

(I-74) 

-A"  /jJwZy;  +sV"°'/3J(u«//3'  Z&*//t  + 4i*fy 

~/V"°'/3cfW/0,/i  / /V^JC  U^//3  / ^ 

S'  M / J~ $<*//3  S-  Q J ( ^ / /t S^)J  a//)  . 

Dy  Equation  (1-43),  the  theorem  of  minimum  potential  energy 
yields 

/ fz/TSu^, 

/ Arf’Cti.  Sty  * - ( Af"?*- sr’f’M  ■ Mfy, 

+ a*3fe«  / - fSu/J  <y/>  = 0 . (1'75a 


Py  application  of  the  divergence  theorem,  Equation  (l-75a) 
becomes 

^ ^ '/ ^ <=V 

- <?■";  / «,  </-*  ^ ^ ^ -/ /Ar 

f(/V'  (l-75b) 

/W-0  / fs  [/!/“/’(*&  *fi,)?(sra'/i+/v~’?3jA< 

- Q^Js*  i)f  dS  fyds  = o , 


where 


^ o//i> 

At  ' 


£(#'-P-  /v"A  t 


(1-76) 
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Since  / u ^ and  /w  are  independent, 

from  Equation  (l-75b)  we  obtain  a set  of  equilibrium 
equations 


O 


(1 


- Q**  = o , (l 

0 

* Q-'U*  + f = o , 


and  boundary  conditions  in  terms  of  stress  resultants 

/rf  % <r  = o , (1 

= o . (1 

/*/"'%  - QPS] )j,  <5w  = o , (1 


( /V' 


t> V'J* 
/ 


/ At*  ' //«&,  = 


-77a) 


-77b) 


-77c) 


-78a) 


-78b) 


-78c) 


(l-78d) 
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The  stress  resultant-displacement  relations  can  be  obtained 
by  using  Equations  (l-50a),  (l-50b) , (l-50c) , (l-55a),  (l-55b), 
(l-71a),  (l-71b),  ( l-73a ) , and  (l-73b) 


/V"  /4>-  /)  ft( +2  J (1=79) 

+ T'*/3  , 

where  T 5 01  ^ is  defined  by  Equation  (l~56b), 


where  T"  “ J is  defined  by  Equation  (l-57b), 

■*  ( “%  **>**%) 

'*%*]  * ■ 


3 


with 


1 

Mr'  =/. 


yx 

2 


ft 


/W 


2f'/iT~(*M  * j)  ^ - 


^ i + ir  Hj/yj  ] V ^ 


(l-8la) 


(l-81b) 

(l-82a) 


(l-82b) 
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//  / - d ^ t (U>/J  f y 2 ^ (1-83) 

— o/ J3 

where  T is  defined  by  Equation  (l-58b), 

/7^=  A ^yJ  pf  6 1-  Atp*  , (1-84) 

— d/3 

where  My  is  defined  by  Equation  (l-58d),  and 

Q*3  = + w)  . (1-85) 

Also,  from  Equation  (l-53a), 

/v“f  = t(^u^  + 

t(  u /yi  + ^w/j  i- j Hpr> %)■?  (1-86) 

-r^/3  + T*°'4  + r«p  t 


and,  from  Equation  (1-76), 


'f  = S- j - t*ms)+  T'^ 
- 7**]  i-x-r^ge*,  / m?  . 


(1-87) 


Py  substituting  Equations  (1-86),  (1-87),  (1=85),  ( 1 -81  a ) 
and  (l-82a)  into  Equations  (l-77a),  (l-77b)  and  (l-77c),  we 
have  the  corresponding  equilibrium  equations  in  terms  of 
the  displacement  variables  u ^ and  w.  The  boundary 

conditions  can  be  obtained  similarly.  The  practical  appli- 
cation of  these  equations  will  be  discussed  in  Chapter  II. 


CHAPTER  II 


PRACTICAL  APPLICATION 

The  above  derivation  is  general  enough  to  apply 
to  any  coordinate  system.  From  a practical  viewpoint 
the  coordiante  system  is  chosen  in  accordance  with  the 
actual  problem  in  order  to  save  time  and  labor.  Here,  we 
consider  an  axisymmetrical  case  of  the  sandwich  shallow 
shell.  In  this  case  it  is  most  desirable  to  choose 
cylindrical  coordinates  r,  0,  and  z.  Also,  due  to  the 
axisymmetrical  nature,  the  displacement  is  independent  of 
the  coordinate  0.  In  the  following,  we  are  going  to  convert 
the  tensorial  components  into  physical  components  and 
the  equations  in  tensorial  form  into  engineering  form. 

We  list  the  non-zero  components  of  the  metric  tensor  and 
Christoffel  symbols  in  a cylindrical  coordinate  system. 

Thus 


<^„  = /, 

= r , 

e"  = /, 

e*2 

/ 

- r*  , 

(2-1 ) 

rA  =-r. 

r2 
/ /2 

=ri  = t . 
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1 . Physical  Components 

A.  Physical  Components  of  the  Elastic  Coefficients  and 
Thermal  Coefficients 


Like  most  problems  in  shell  theory  the  elastic 
coefficients  are  obtained  from  a generalized  plane-stress 
assumption. 


A 


////■ 


For  Isotropic  Material 


/-  >>* 


A 


//22 


A 


/2/2 


. £ A __ 

'/-}>*  A 

.,2/2/  . 2H2 

/)=A=A 


£ 


2(/ti»  / 


A 


2222 


/- 


(2-2) 


JD 


// 


A 


22 


/-  J 
2 o< 
/-/> 


(2-3) 


where  E is  Young's  modulus  and  9 is  Poisson's  ratio. 


— //// 

A 


b.  For  Orthotropic  Material  (30) 
£r  _ 2222  £». 


'-A*,  A* 


_ 2222 

A 


At 


-//22  -22//  __  £» 

A A A-fi- 


j/2/2 


£ 


A 


-/y/y 


2 , 


A - 2 A 

A — 22  a?, 


(2-4) 
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n 7/  - - ^ J<ba(&) 

V ~ /-  ^ >A 

A f ^Sl  ^ "^Sl0*/  ^ 


(2-5 ) 


/ - y y 

' -V-SSl  -*»/■ 

B.  Physical  Componets  of  the  Derivatives 

For  the  first  covariant  derivatives  it  has  been 
shown  in  (27)  that 

o'  </* 


3 -*X> 


^ /v  * > ' s> 

From  Equation  (2-6),  we  have 


A'  4 

= A/", 

t 

y 

aa 

- /^, 

-h 

/;> 

/2 

} 

/ V 

= rJ, 

*/;/ 

A 

— p r* ' 

/V 

-?/ 

j 

= ^v 

/ 

/ V**  . 

= y/7/ 

A" 

/<? 

(2-6) 


(2-7) 


For  the  second  covariant  derivatives  we  obtain 


o//? 


J / 'CO  J 


-C#*  • (2'8) 
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Applying 

formulae 

(2-6), 

(2-7) 

and  (2-8) , we  have 

fa 

\ 

\ * 

\ 

* 

II 

/ ' 

fa 

=rAfa  + nl 

/ / 

Ws{ 

(2-9) 

fa 

= 2 r 2 

fa 

- rj;  fa  , 

= fa* 

cyQ 

' fa  fay 

(2-10) 

From  Equation  (2=10),  we  have 

fa  = . 

fa>2  — - f**  fa  > (2-11) 

fa/  = M//  . 

Similarly  Equation  (2-11)  can  also  be  applied  to  Z . 

Also,  we  have 


Hence , 

^ 2/ ^ ~ 


l-y 


(2-12) 


(2-13) 


C.  Physical  Components  of  the  Stress  Resultants 


The  relation  between  physical  components  and 
tensorial  components  is  expressed  as  follows 

p — fP  p*/3 

'(*{*>  r ' (2-14) 

p^>  ~ JP  PP  . 

From  Equation  (2-14),  we  convert  all  the  tensorial  com- 
ponents into  the  physical  components,  thus 

' /y"  ~ A/f  , /V" 

, Af“^/Wr,  (2-15) 

A/3*  ~ • 

2.  Equilibrium  Equations  and  Boundary  Conditions  in. 
Engineering  Form 

We  use  Equations  (2=7),  (2-9)  and  (2=15)  to 

obtain 


<AA/r 
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A.  For  Method  A 

T 'sing  Equation  (2-16),  we  can  have,  hy  introducing 
two  new  variables  X = i ^ + u"°<  ^ and  ^(u^  ” ) , 

the  equilibrium  Equation  (1-52)  written  in  terms  of  the 
engineering  form  of  the  stress  resultants.  Thus 


r 


c/  /Vr 

y- 

- ^ 

— o 

c//Wr 

& 

A 

- 

cVr 

-?r 


'/'A'  ^ i- 


(2-17a) 


where 


4 = /*/,.  ^ 

3*  = ^ ^ . 


(2-17b) 


a . Equilibrium  Equations  in  Terms  of  Physical 
Displacements 

The  stress  resultants  can  be  expressed  in  terms  of 
the  engineering  displacements  as  follows 


(681-3) 


( J8T-3) 

(981-3) 

(P8I-3) 

(381-3) 

(qex-3) 

( e8I-3) 


V ^ [('4?  f +%>)-r  + (*0  / + 


S/jL. 


< *00  +[Wf+  *)?$  * 

JP  li'Tt*!' -/)3  rj,  , /a  1 tr 

+ XSrV^-  * — 


?1? 


t<  frpJ  ?J7 


' [ffpri  + ij?  / ^ J'£ 


7 


( 74  [(^/*  f 7 */?  )stfc  ^ ^0  f 


, sr  k^^-09  v 
tfrJ  “ J ~*7 


-/  , rSP  ] (7  -09 

O AfrpJ  f)7~ 


*(7  +fr&f+  '*)-?+(&§ 


c 

s/*l*fr  7)9  , r y , ✓ D,« 

WJ  V rf*3  l wr  + zrp  <* J-ji 


(7-0_, 

7 - 


*/r  -/ 


Aj  +(##+/#)*+&)  rj—j^ 
*7[j  +/&)?< 

•'_/ .Td-'*  *■*£■*£ -L  f'm  , '/>]'*•*■-/ 

M 1 J ' MS’ ZP  JmF/t  + ?£r  J —pj- + 

S £*  1 rt  -/ 


/ ^ /<  ?'  ^7/  , * A//?  ?/?  + Jj f. 


)7Z 
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where 


Tr  “ T/  + T,"  + T,  , 

(2-18h) 

' » * / --m  - //  ~ 

/ tsl  ~ n&.  A 71  /■&.  j 

( 2-18i ) 

(MT)r  = ( A,)),  r +(Mr),  , 

(2-13 j ) 

SA,  — i^r)»  r (^rJ«  f far)#  ' 

( 2-1 8k ) 

= (A7r)r  +t(T/-T/)  , 

(2-181) 

(A)*  = (At)„  + i( A- Tj.)  . 

(2-18m) 

Py  substituting  Equations  (2-18)  into  Equations  (2-17)  we 
have  the  equilibrium  equations  in  terms  of  the  engineering 
displacements 


y / *2fr  . c/JV  dT  dJd  , dd  d*d  ),  dfr 

• I dr*  dr*  dr  dr  dr*  r dr  dr*  s'  dr 

+*;f. '*(&)%  5?3tJ  -4  £ + a:  (/  *7' 

i-Tr  -71)=  O , 


dr 


(2-19a) 


(2-19b) 
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<=/  Jr  <=/Z  J.  { j 


c Jr  * ‘Jr 


Jr 


^ * ( rf  fs  ^ ^ ' '/o 


< zh 


(r&  + r&)  +fa  (0  / 00  ^(0/J 

+ + 4°  71  Jf '0  ? 0J  / 4'afr  / 

0,(r^ 


a'*W 


c=/r*  ^ ^ 0/'J  J * sjrj  ( 4 


7 c//V  _/  ^V, 

I 0 ^ /■  ^/v 


7 4o  7r  ( r) 


where 


^/  “ y 


' ^ £,t(/'  J&r)  J / 


^ - 


_ f 2£t 


/-  »* 


t 


4*  £ 


4. 


(2-19c) 


(2—1 9d ) 


/ -4s 


(2-19e) 


(2-19f ) 


£^_t  t* 


£tt 

/-y 


v- 


£+  t * ^ £ tt 

6('-~*rAr)  /-»*  (2-19g) 

£-t  t Sj,  £*  } 


2 &t 

£ tt  £,  t* 

f-»*  f £(''>UV 


(2-19h) 


£+tt* 

/£( * ~ 

£tt  ^ £r  t * 

''»*  ^ 4r) 


(2—1 9i ) 


As  ( t + t) 
2 


(2-19j) 


£r  ? ) / 

/ " Ks  ^9/  J /0  ’ 


zE  tit 
/-  ti2 


(2-19k ) 
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+ ,)y 

1 ^ /0 


Jf'  = ' 

^ 2£t  + £r  t 


£tt* 

/2^W 


£ t* 

£(/-  »V 


£ t ( t ft) 

/-  F 


£t3 

<o  ( / - v V 


y- 


£ t(  ttt) 

/-y* 


F-r(r>=  * 


// 

c**f#rJr 

Jy  * 

cy?* 

■*l 

r , d*(Wr)r 
' Jr* 

c/C'VrJa 
cJr  * 

* 4^]. 

(2-191) 
(2-1 9m) 


(2-19n) 


(2-19o) 


(2-19f) 


(2-19q) 


(2-19r) 


b-  Boundary  Conditions  in  Terms  of  Physical 
Displacements 
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^he  boundary  conditions  can  be  expressed  by 


A/r  J Jr  ~ O , 

(2-20a) 

= O , 

(2-20b) 

( r/V,a,(»"z)  + + /*,-/%,  t(r ^ 

+ /f,~  At*)  + t '<?,  //V  = 

(2-20c) 

(2-20d) 


Ry  using  Equations  (2-18a),  (2°18b),  (2-18c),  (2-18d), 
(2-18e),  (2—1 3f ) and  (2-18g),  we  obtain  the  boundary  con- 
ditions in  terms  of  engineering  displacements 


J J J/y 
Jr  jr 
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/ Jm/  \2 
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(2-21a ) 
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(2-21b) 
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dr 


d / c/ m/ 
dr  dr 


z I dr/  J 
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c/ ur 


oY '*fr)r 


1-  4.  r 5? +lr  + tor)r  -MU 
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d JV 
>/  M/-2 


a'  ~L  djd 
T ^ A ma 


where 
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d 1 1 


/-  r* 


dr  t* 


^ ~ /-i/'  ' 


^r_M»r_  ^ 


(2-21c) 


(2-21d) 


(2-22a) 
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(2-22b) 
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^7 

_ zr,.  1 1 

10 ( 7 ” 

) 

(2-22c) 

^/d 

- -A  tt 

/ 

(2-22d) 

9 

i 

- 2^(  t j 

; 

(2-22e) 

& ( t / Tf 

t 

j 

(2-22f ) 

_ £t* 

t t* 

(2-22g) 

&(/->*) 

s»tJ 

£,  ii,t  t* 

(2-22h) 

B.  For  Method  B 

Substituting  Equation  (2-16)  into  the  equilibrium 
Equations  (1-77),  we  have  in  terms  of  the  stress  resultants 
in  engineering  form 

K ^ SK  ~ ^ - ° ' (2-23a ) 

XI  XI 

r + Mr r<?,  = ° ' 


( 2— 23b ) 
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SVr  ( dr* 


i- 


dJd 

dr2 


) + 


r l dr 


+ Mr)  7 r dr  (A/;  +M*)  ~ r-Jp{A&  rA^'j  / fi£fr<p,  ) 
+ ? = ° • 


(2-23c) 


boundary  conditions  become 

/A,  <f  U,  ■=  O , 

(2-24a) 

r\ 

A/r  je,  =o  y 

(2-24b) 

[ r (i?  * +£('/%  * r?p  + dd'd-Ati 

- a£)  - Qr  U»  = o - 

(2-24c) 

( a*;  / At;j  s(0j  = o . 

(2-24d) 

a-  Equilibrium  Equations  in  Terms  of  the 
Physical  Displacements 


Py  a similar  procedure,  we  are  able  to  express  the 
stress  resultants  and  moments  in  terms  of  engineering  dis- 
placements as 


d£  t r ^ dddjr 
/-  j,*  l dr  r dr  dr 


£ f dS  , d-  / c/dA 

l d'r  / dr  dr  r * ( drJ 
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rr 


(2-25a) 
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where 


7T  = T/+  TrW  T,  , 

7^  = £'  * £"  r&  . 

^Vc// 
c/r  c/r  c/r 

W ^ t‘'rdr  'V dr'J  r r I '*■ 


Ai,  = TTA'j  A l dr  * dr  dr  ^/(drjJ^  7 j 


2>£t 

v»  ~ /-> 
/-  rsa 
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From  Equation  (1=87), 

^ 7T  / f - d&r  ,*  £^*/  ) / - &r  , * / <S/V,7 

A =77-711  u 77'(  TP  T ^ ~ rsr^ 
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'ZO-JinM  I e/t 


/}  [ ./.TcJ&r  i/  //  ^ / gV,/ 
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(■ V*  ^ + 7j  + (/$■, 


7~'&. 


From  Equation  (1-85) 


Qr  — <7,  t ( <9,  -t  yy J , 


and  from  Equations  (l-81a)  and  (l-82a), 

, „ JT  f*  f — ' 

Mr  -f  A/,  = 2(/-j)‘j[  ^ dr  J 


-?rf9]f^rJ,+Cvy  , 


(2-25b) 
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(2-25d) 


2-25e ) 


(2-26a) 
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-*J*  [t 


ft 


Jr* 


-A 


(2-26b) 


Substituting  Equations  (2-25)  and  (2-26)  into  Equations  (2-23) 
we  obtain  the  equilibrium  equations  in  terms  of  the  engineer- 
ing displacements  as 


cJ  S'* 
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d J, 
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, J J JS 
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(2-27d) 
(2-27e) 
(2-27f ) 


(2=27g) 
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A?  - 2(t-»2)  * /2( /- ^r9?ir) 
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E t t* 

2(/-j>V 


y^7  ~ <?,  t 

* *£tL—  + 

~ s(/-y*j 


£ -t  ^ 


2 ( / - y*J 


(2-27h) 

(2-27i) 
(2-27 j ) 

(2-27k) 


b.  Boundary  Conditions  In  Terms  of  the  Physical 
Displacements 

When  we  use  Equations  (2-25)  and  (2-26),  the  bound- 
ary conditions,  Equations  (2-24),  can  also  be  expressed  in 
terms  of  the  physical  displacements  as 


r r dnf  . d(/ d?  . m ( dhi  iJ  / j.  / a _ /t  ) ^£- 
/ A, I dr  + dFa?  * d Id/)  ) r 

+ T,  } = ° , 
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t/e,  A y.  J-  cM 


(2-28a) 


^ f ?%?)+  Jr'r'#,  ^ 


, 4?  Er'  '4  t ^ 
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T,  }($  + &)+ +>?)-&,-*)(£ 

/ jr  0?)  + (At'A  + (/Or  J + 


(2-28b) 
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(^0  ~^e)(  dr*  / dr'j  ~ ^7  ( / 7y * 


-(A/tJ?]}  cfh/=o 

(2-28c ) 

U(i?  +?Q-(4-4)fcZ  -/^^^//^(2.28d) 

where 

/)  _ ^ v ^r  t 

(2=28e) 

0 

' r ZO-^J 

The  conditions  at  the  center, 

r = 0,  due  to  the 

axisymmetrical  nature  of  the  problem, 

are  as  follows 

Ur  — O , 

(2-29a ) 

Or  — O . 

(2-29b) 

dW 

dr  ~ , 

(2-29c ) 

and 

w is  finite  at  r = 0 . 

(2=29d) 

The  boundary  conditions,  Equations  (2=26),  together  with 
the  conditions  at  the  center,  Equations  (2=29),  and  the 
equilibrium  Equations  (2=27)  form  a complete  specification 
of  the  engineering  problem. 


CHAPTER  III 


NUMERICAL  APPLICATIONS 


In  this  chapter  an  axisymmetrical  sandwich  shallow 
shell  of  paraboloid  geometry,  ^ - r^  ^ _ T^y  is  consid- 

2R0  ro 

ered  as  an  example  for  a numerical  solution.  The  structure 
is  under  uniformly  distributed  mechanical  loading  as  shown 
in  Figure  4.  The  example  does  not  include  thermal  effects. 
However,  the  method  employed  here  is  valid  no  matter  whether 
the  thermal  effect  is  considered  or  not  when  the  temperature 
distribution  is  assumed  as  a known  function.  Then  one  simply 
adds  some  terms  due  to  the  temperature  distribution  to  the 
equilibrium  equations,  as  well  as  to  the  boundary  conditions. 
Furthermore,  we  consider  the  edge  of  the  shell  as  clamped. 

The  paraboloid  geometrical  configuration  is  a good  approxi- 
mation to  a shallow  spherical  case. 


1 . Dimensionless  Parameters 

The  equilibrium  equations  and  the  boundary  condi- 
tions obtained  in  the  previous  chapters  are  in  what  is 
considered  a "dimensional”  form.  Each  of  the  elastic  con- 
stants, E.  p , Er , E@,  PrQ,  geometrical  properties  of  the 
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structure,  t,  t,  r0,  R0,  and  the  space  variable  r have  the 
appropriate  physical  units.  It  is  customarily  desirable  to 
transform  the  differential  equations  as  well  as  the  boundary 
conditions  into  a "dimensionless"  form.  By  solving  the 
"dimensionless"  form,  one  has  the  solution  for  a number 
of  problems  for  a certain  range  of  values  of  the  parameters, 
and  consequently  a greater  amount  of  information  and  more 
general  applicability  will  be  reached. 

Let  us  introduce  a set  of  dimensionless  parameters 


f = 

~k  < / . 

(3-la) 

IV = 

h/ 

h j 

(3-lb) 

u = 

( 3— 1 c ) 

6=- 

/q 

IT  0r  ■ 

(3-ld) 

K ~ 

rV 

(3-le) 

2 /?./)  . 

A.  Equilibrium  Equations  in  Dimensionless  Form 

By  performing  the  transformation  and  dropping  the 
temperature  terms,  Equation  (2-27)  becomes 
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-e<f  + Steffi- rtf ]=o, 

c>  / „ 2(&  . c/& , _ n .&.  _ a />,  c/JH?  , VfJV 
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~ - #7  (?0  + ?if)  = ° ' 


( 3-2a ) 


(3-2b) 
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’ = _ >,-?  <2 W y.  -t- 

vy  y cjy  ^ v/y 
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where 
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(3-3a ) 

II 

(3- 3b) 
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2/hV 
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(3-3c ) 
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4 : 

ll 

( 3— 3d ) 

4 = 

3S 

3g 

( 3-3e ) 

4 = 

^ ' 

(3-3f) 

3?  - 

3?K? 
38  ' 

(3-3g) 

f = 

?rS 

( 3-3h ) 

/)s  /) 

P.  Boundary  Conditions  in  Dimensionless  Form 

For  a clamped  edge,  the  boundary  conditions,  Equa- 
tions (2-28a),  (2-28b)  and  (2-28c),  are  automatically 
satisfied.  Equation  (2-28d)  becomes,  at  1, 

36  _ . q'ty  s 

& d/y  ( / r y J - O,  (3-4a) 

while 

e = a=-^y-o  y=o,  (3.4b) 


6 — u = tv  = o 


(3-4c) 
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2.  Numerical  Solution 

It  is  obviously  very  difficult  to  solve  exactly 
three  coupled  nonlinear  differential  equations.  Accordingly, 
a certain  approximate  solution  has  been  sought.  A perturba- 
tion method  with  the  center  deflection  in  dimensionless  form 
as  the  perturbing  parameter  is  used  to  reduce  the  nonlinear 
differential  equations  into  several  sets  of  linear  differ- 
ential equations  with  variable  coefficients.  In  solving 
each  set  of  linear  differential  equations,  the  power-series 
expansion  method  is  employed.  A general  recurrence  formula 
is  thus  obtained  from  the  differential  equations  and  all 
the  coefficients  for  the  series  can  be  determined  by  imposing 
the  boundary  conditions. 

For  a practical  reason  a finite  number  of  terms 
of  the  power-series  is  assumed  as  an  approximate  solution. 
Fy  using  the  recurrence  formulae  which  are  obtained  from 
the  differential  equations  and  the  boundary  conditions,  one 
can  have  a set  of  linear  algebraic  equations  with  undefined 
coefficients  as  unknowns.  In  solving  these  algebraic 
equations,  a digital  computer  program  of  the  matrix  in- 
version written  in  Fortran  II  has  been  used.  The  computed 
results  are  plotted  in  Figure  5. 

The  problem  is  then  reduced  to  an  ordinary  shallow 
shell  by  taking  the  thickness  of  the  core  equal  to  zero; 
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i.e.,  t = 0.  By  doing  so,  one  merely  has  two  facings  bonded 
together  working  as  an  ordinary  shell  structure  of  thickness 
h = 2t.  The  resulting  curves  are  compared  with  known 
solutions  (23,  2A)  as  shown  in  Figure  6. 

A.  Sandwich  Shell  Structure 

For  an  isotropic  material  we  have  B4  = 85, 


rl  - b2*  r1  - ~ b2  - (2B3  - 4Bj)k2  and  B4  = B3  = B4  + B$. 
heref ore , Equations  (3-2)  can  be  rewritten  in  the  form 
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f 8,  k[ 70--  0{0)*  tuJ-B,ku 

0(0f]+ 

■t  B,K  y } 37  / 0,  [?cf  - 7*0  7 (07 

W - (70 - f 0 ' r^fm 


7 — y y.  0-7  c/0  0 1 — 

< *?*  ~ 7 0?  + yt)  + f 7 = o . 


(3-5b) 


(3-5c) 
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Let  the  functions  W,  u,  0 and  q be  expanded 

in  the  follow- 

ing  form 

IV  — IV/  7>  w.  y wf  + • • • , 

(3-6a ) 

U — U/y)  iv.  + UJfiw*  / • • • v 

(3-6b) 

6 - 0,(y)  K -f  Q(y)  tv.*  / - - - , 

(3-6c) 

f = //  *•  + 7>"*2  + 

( 3-6d ) 

Substituting  Equation  (3-6)  into  Equation  (3-5)  and  equating 
similar  powers  of  w0  result  in  a sequence  of  sets  of 
simultaneous  equations  for  Wn(^  , Unf and  qn. 
Considering  on  expansion  of  two  terms,  we  have  two  sets  of 
simultaneous  equations.  From  the  first  set, 


B,[  7 


4 ? 


c/V 


¥-(r* 


*7 


— o 


a r + /A  - v.1  - 3 f # / — 

- faf)  - 4 (T*.  + ?0->  = ° - 
4 (r#  - f/  - u)-(r0**0-  - 


( 3— 7a ) 


(3-7b) 
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_l  c/e,  , e, 


7 <*7 


+ %)  + ?.  7 = 


= O 


'he  corresponding  boundary  conditions  are 


d,  = U,=  =o 


IV,  = / 


f 


= o 


and 


0,  = U,  = o 

U <70,  _ 

6 7y  ' ' c^/{  <7  y*  ‘ f 7y . 


% si- -('-%)( =° 


"he  second  set  of  differential  equations  are 


b[  - Ul  _ /#*</% 

a'l  r *?*  + ? (?  -7^- 


(r&  +?&]- 

+ £ r(/f)(0)  / £@?f=  o , 


4 (70-  '0-f-J-4(70  + 0 
- f~  - 07(?q.  / = ° - 


{Zk(70-  f0-'W-£jKi/J0  / 

+ f <4-  *>■>* (0-  - ?0.)  + 


(3-7c ) 


(3-7d) 


(3-7e) 


(3-8a) 


(3-8b) 


■ 7(  <*? 


,jA 


i-  0 ^ % -L  cSQ  &L  , ^ 

j("^r  *?*  ~?  *7 ° 


(3-8c) 


he  corresponding  boundary  condition 

u.  = e<  = sf  = hi  = a 


s are 


*#/ 


7 = ° * 


(3-8d) 


and 


u*  = e*=H£  = o 


7*  7 ^7  ■ 

y - / . 

Let  us  solve  the  first  set  of  differential  equations. 
Equation  (3-7a),  we  have 


(3-8e) 


From 


7 ~ w + <?,  ? . 


(3-9) 


vhere  c\  is  an  integration  constant.  Py  substituting  Equa- 
ion  ( -9)  into  Equation  (3-7c)  and  performing  an  integration 


64 


with  respect  to  ^ once,  we  have 

_ jg,  __  _/£¥ 

^ ( 7 <s/p*  f cJy  p J ~ ? <z/?*  * ^ pi  ~ ? c/f 

- 4s  ?*+  *r.7‘  ■ 


-Q-  ^ 


(3-10) 


Substituting  Equation  (3-10)  into  Equation  (3-8b),  and 
solving  for  G^ , we  obtain 


a = 


7 JJx  _ ~ , J-  <^K\ 

&7  l ^ QJl  <=/?*  ? cJp*  ?*  7 ' 


- ^ 


3*  £-1 


L _ B.  0*  #4  X 

£7  7 3,  3j  3y  ?• 

_4/'  - cJrt 

ZSj  B7  7 " ^ 7 ■ 


£0* 

c'7 


(3-11) 


Let  the  function  be  expressed  in  a power  series 


00 


(3-12) 


By  substituting  Equation  (3-12)  into  Equation  (3-11)  we 
obtain 


« =-i 


4_^l 


4*7  7 


r / 4 

^ l ( Bj  By 


- 


- -&.&*  c 


~ 2 


‘Jr  +£,[(& 


§*- ) £('/>  t o(/>  ftfot  ~ -■*£- 


- 2/'///^^  ] 


7 


*?/>  y-/ 


( 3=1 3a ) 
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By  applying  the  boundary  condition  (3-7d),  we  have  02 
and 


9,=  ^y 


where 


A 


*>*7  A 


J Of".  - &! 


7 2 &,  3j  &7 


_ *if+  c 

**3*?  ' 


_ ■? 


&3  e7 


and 


- _ o /*  /•/ ) o/ 

■f  = t.  2.  3.  ••  - 

Also,  from  Equation  (3-9),  we  obtain 

« = £/v.  ^". 


where 


7 


and 


f 3,  J > 


c 


= /.  ^ 3 \ 


0, 


3-13b) 


-14a) 


3-14b) 


3-15) 


3-16a) 


3-16b) 
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Substituting  Equations  (3-12)  and  (3-13a)  into  Equation 
(3-10)  and  equating  the  coefficients  of  the  same  order 
terms  we  obtain  the  recurrence  formula  for  the  oCs  as 


&.  0* 
2 0, 


J2(  / + + //S2(- 


_ 


V 


- 3 gJJ^  - 2 M = o 


y — 6 


7 


7' 


and 


/6('r  - 4 fy  * 

= o . 


-fA3.  cv 

e, 


Q0*  _,y , 

S7 

^ - 7.  3,  4. 


(3-17) 


(3-18) 


py  this  recurrence  formula  (3-18)  it  can  be  proven  that  the 
power  series  is  convergent  (Appendix  A).  We  have  only  four 
unknowns  d 2,  ^4,  and  q^  to  be  solved  from  four  boun- 
dary conditions,  Equation  (3-7e). 

As  a rough  approximation,  let  us  assume 


2 


/ 

7 ^ ^7  , 

(3=19a ) 

(3=19b) 

l7r  + ■ 

(3-19c) 
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Py  using  the  recurrence  formulae  we  have  a set  of  algebraic 
equations 


~jf  ~ ^ 4?  °(z  - 32  ( / y 3jJc2^  + / / S2  4 e£  / 

1-  A c - ^ 3, 

* ' ~ 2 0, 

( 3~20a ) 

°<a  ^ o<V  y - ~ / , 

(3-20b) 

&3};  y-  *4^  ^ 

(3~20c) 

>%  y £ / 74  £ -(/-£t)[j(u 

^ /jVo^  y ^ruTy  ^>0/  ]=  0 , 

(3-20d) 

f3'  '3  "A  *^=0, 

(3=20e ) 

fa  ~ y , 

( 3=20f ) 

N 

0 

1! 

t 

4' 

Nh 

* 

(3~20g) 

fJ?  f * **  ^ O , 

(3=20h ) 

/7  _ J=L  — 

/ " ^ ' ^47  / 

> - 4jlA>  , jl  ^ 

4 V ' ^"4  4 ^ ^ #jB7  + 

(3-20i) 

/ i’  c/,  = - 4? 

7 

(3-20j) 

= 


% 74  % + 4 ^ 4 


( 3=20k ) 
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ys  + z ^ = °, 


(3-201) 


( 3=20m) 


where 


4*  _ 


_ & #4 


s,  b7  > 


(3=20n) 


These  algebraic  equations  (3-20)  are  solved  by  the  first  part 
of  a digital  computer  program  in  Fortran  II  (Appendix  B). 

The  results  are  to  be  used  for  the  second  set  of  the  diff- 
erential equations. 

The  same  procedure  will  be  applied  in  solving  the 
second  set  of  the  differential  equations  in  which  the  non- 
linear terms  of  the  W]_  and  Uj_  are  treated  as  known  functions 
of  ^ . From  Equation  (3-8a)  we  have 


(3-21a ) 


where 


. ( 3=21b) 
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Py  inte9ration  of  Equation  (3-21a)  with  respect  to  ^ , 
obtain 


we 


*7 y ^ - £^f?'  / 


^ <5  / , 


f 


(3-22) 


where  is  an  integration  constant. 

Substituting  Equation  (3=22)  into  (3-8c)  and  integrating 
once  with  respect  to  ^ , we  obtain 


+ f=i(70-  + ^p. _fg}_ 

~ ■ tj 

(3-23) 


where  C2  is  an  integration  constant, 
and 

^>o 


# r =j4F  -^/W ■/?  =ix  ^ "'VZ'  £ 7***-/' 

-(pts-OBj  -JAs(f-o[ 

2(23  + 00, 


(3-24) 


Substituting  Equation  (3-23)  into  Equation  (3=8b)  and  solving 
for  ©2.  we  have 
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£ 


*3 


_2.  _Z  ) _ 

* ? <=2 fr*  ?*  o'p~s 


f <J?J 

&*  , 


' 4%  4 rf?**^?* 


t &A 

A 


t£+  wwwr-m  *?- 


_ 


— 7 „ A0£  _ 

2Ar3j  1*7  c=S? 


fs 
As  7 


Let  W2  be  expressed  in  a power=series 


po 

2: 

/=/ 


Wx  = £ ^7  T*** 


Substituting  Equation  (3=26)  into  Equation  (3=25), 


* %f  fc,  ^ 


where 


J2  &,  — A*  / 

% ~ gj  4 (4£,  ~ AAj)A/3j 


- A0, 


-=z-^=—  c _ _r_  9 - ^ of 

A0j07  ' ^44  T*  * 


-/  &j  c=^>  ^ ^ ^ J 2 f>  (*£,  - 

- £3  ) A 

-0SV-0JA,  -A&]  ^(^-s^/J^s-2}}  ~ 


(3-25) 


(3-26) 


(3=27a) 


( 3=27b) 
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- A.  , B,  &+  Q A _ 

^ B,3jB7K  ?(?-/J(  * < r 


- 2S-0 + 


°A/  -jy 


<=^ 


(3-27c) 


for  p = 2,3,4 

and  C2  r 0,  by  the  boundary  condition  (3-8d). 

Substituting  Equation  (3-26)  into  Equation  (3-22)  and 
solving  for  U2’ 

64  = 5 /?  (3-28a) 

^ /-?/*•*/  / 

where 

/3  = , (3-28b) 

and 

/*_£_  , Bj  7 - _ / s(t~s+Q 

<23 -i- fj  + °(z(f-s+/J  , (3. 28c ) 

for  p : 1,2,3,  ... 

Finally  by  imposing  Equations  (3-26),  (3-27)  and  (3-28) 
on  Equation  (2-23)  and  equating  the  coefficients  for  the 
same  order  terms,  we  have  the  recurrence  formula  for  c(  9 s. 
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//S2 1,  ^ y ^ 

= *§t£[g(££) -gj)/}^  + -j{*4- 

- *(**-*)*fi~,  - v4a: 


y £sA/,[(j6. 


S=2 


s y-f, 


2(?-*+*)s%-(r+»0t]o<,(/.J[„J/3MS_ 

- f (J-vfet 2 f -SS  fSjB,  -JBj  Ja^ 


/ 


^2S-2 


'*?(' /&-S-S2) 

_/r 
/ 


- ?£,]°<s(Jt_Sf0fis_r 

- (yi-s*o(s- o[?( Jfi-BS-t- VB,  -J&j '.]°ty.jJ 


o< 


2S~2 


}%z£  sf?  -*">/*(*?  * 


+%] 


IS^Cf-s.o 


(3-29) 
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-2S-08,*  B3]ot;t 


(/>- s) 


(3-30) 


for  p = 2,3,4,  . . . 

Here  we  have  four  unknowns  <*4,  c1  and  q2  to  be 

determined  from  the  four  boundary  conditions,  Equation 
(3-8e ) . 

We  set  approximately 


^ ^ ^ ^ ; (3- 31a) 

<2?  = ?7+ <3-31b> 

^ =/?r  *fi7J+/?ir7si-/*777i-fe?9i-fi,y'-  (3-3ic ) 

correspondingly,  we  have  a second  set  cf  algebraic  equations 
to  determine  these  coefficients 


dr  - ^4=>i  - ‘4+ §-5 

--  ~L[ 3(20,  - 4/*5  - 4,)<=<,fo  - 

- 3(*Z  - BjJo/Z] - B,y, P ^ - 


2 

'<*T 


(3-32a) 
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(3-32b) 


4 £ ^ ^4  °£  ~ /**  l,^  + ^ ^ 

= - *+* 

B7 


[*(■**.- sj)fi<x*  *■ 


3(-*e, -&,)<=£] l + 


3,B« 


& . 


^ % 


(3-32c) 


%[>,  1-3%  +S£  t7>-7  f 9%  +■//%]-< '/-Bt)[ 

2(/fH)37  / £(S fJJ37&  /■ 2f 7fJ>)^  i-/0( 7f3jjJ- 0,  (3-32d) 


ft  + fs  * f3?  * f?  7-  f„  = O , 

&)<**  = , 


/j-  ~(j  f60,)  °<+  — J£,/<  (**!  * , 

f^y  ~ ( ^ ^ 5 0,  ) ^ — 


* 50,  / <=*6  ~ #fAr  i l ^ j ,-  .^~6 


+ #>sU  * 


(3® 32e ) 
(3-32f ) 

( 3® 32g ) 


’Y#’1  #></,  (3-32h) 


c. 


— o 


(3-321) 
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^ 4 ^ / Zp.^~B7  ? <2  3s07  f 2 ^ 


°<2 


- - 


4^7 


^ ^ 4 ^ £ + >7  + >9  + t,  = ° > 


k _ 


4 - 


= - _4^ 


4 ^7 


' f&(^B<  Q)^,  t4  y [^-(f#, + 

jj  /• 


> 4f< 

^§4* 


r ^ * -§-j^< 


* j 


y7  /«  i ^ ^<9=^ 

= - ^Xr 


— ^ X"  r 

^r^l24(33'-BJ^6/?j  +/£(34-%)°<tfr  + 

-/■£(<?#,-  £j)^^7  -4(/o3f-J3J)<x£-6(/o0t 

3 ^ (JB,^  §-)(?- -f 

3 *~o/+. )] , 


(3-32j) 
(3- 32k) 


(3-321) 


(3-32m) 


92/6  ^ ^-4  ^ ~ /4£(/ J- t ~~  c=23 

~ “j^y  /^^?4  - 4-^=^  /•  /efays,  -JS,)^  + 
+ J(2#- - 2(22#-  jf  - 
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- f 4 (XJi,  - fa)/<fa 

~ V ~ ^Ts^<  ^ b>  * y i 

3X400  4 ^ “ ^4°S  " 3X4-0+  %)<*&  y-  fjjf  cfa 

- 34  x 3^  r,  _, 

/ t ■ ^ "4^  ( '*  2 <=& fa  ) / -^4  “ 

- fa)c^fa7  - £ (/o B,  ~ 3 3j ) <=y  - / * ( X #,-33j * 

<**]-  *(*4-%)X/i,**s  - £[+(**,  ~J^fa 
/ ~ ~ XX,  -3Q  J- 

*^3/  +**)-£*•  * (6*,  £ht(*4  + 4M 


4 


4*+  . 


oS 

■*  J 


- ^ ^ ^ 

A ' ^ y‘7Z&)Q<t  = ~2^X&,  (?#, 


* j 


f3"  ~ ( Z~  ^ '/*#,')  fa0  = -^XTFb,  (/0%  * &j)fa 


2 

6 / 


y ^ _^£_  — , — A/r 


s^j)^fa^  4 4( /6  3,  - J-3j)*yfa7  - 36(43,- 


] 4 ~-fa „/^z^> 

^ 6-/  i 


-y 


* 6 ' 


3-32n) 


(3-32o) 


3-32p ) 
3-32q ) 


3-32r ) 
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(3- 32s) 


These  algebraic  equations  are  solved  by  the  second  part  of 
the  same  digital  computer  program  in  Fortran  II  (Appendix 
B).  The  results  are  listed  in  Tables  1 through  4 for 
various  values  of  k which  is  the  index  of  the  shallowness 
of  the  structure.  The  dimensionless  loading  versus  the 
dimensionless  central  deflection  curves  are  shown  in 
Figure  5. 

B.  Ordinary  Shallow  Spherical  Shell  (as  a limiting  case) 

If  we  set  the  thickness  of  the  core,  t,  equal  to 
zero  the  case  is  reduced  to  that  of  an  ordinary  shallow 
spherical  shell.  Thus  B4  r B5  = B6  = B7  = 0,  B3  = B4  = 0. 
The  number  of  the  differential  equations  of  each  set  is 
reduced  to  two.  For  the  first  set  of  differential  equa- 
tions, we  have 


(3-33a) 
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*<[ ~ +o,]-(. ^0 


where 


& = ( 


and  the  boundary  conditions 


W ~ /y  U/—  dy  ~ ° ??  = 


y = *v,= 


- dM4  __ 


7 


<*/-  y = / 


The  second  set  of  differential  equations  become 
d L/j  d 0. ? __  ^ d dHj, 


*(r*& 


_ f -y.2  <d  ) 7 

I ^ 7 <dy  JJ 


<dy*  d y y l { dy 

■+  j.  / d d d/Vj 

■?/<  I d y d yi  T ^ ( cd  y / ~ f dy  ~~  ° / 

' + Ut  - T‘0  rmy]-(?0  . 


d y ' 

-/  ^ °LJM  _ J_  d*A&  _/  c/ttA  , ~ 

dy3  f dy2  ?*  ) + f*  ? - ° 


(3-33b) 


(3-34a ) 


(3 -34b ) 


(3-35a) 


(3-35b) 
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and  the  boundary  conditions 

are 

U2-  dJf  - )V2  - o , 

«=>/ 

7=  °. 

U2-  -M-o, 

7*' 

From  equation  (3-33a),  solving  for  U^,  we  have 


<-c,p  . 

Substituting  into  Equation  (3-33b),  we  obtain 
W S-  c,  7 ( ? 7-  £ - y ^^2  ^ 


B,  Bj 


+ r + ?>r  = ° • 


Let 


=e  . 


^ 2i,-4 

~X~?7=2^ 

7^ 


(3-36a ) 
(3-36b) 

(3-37) 

(3-38) 

(3-39) 

(3-40) 
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oubst ituti ng  into  Equation  (3-38),  we  have  the  recurrence 
relation 


^ * B,ct  + = 


for  p = 1 , 2,  3 


Substituting  Equation  (3-39)  into  Equation  (3-37) 

•*JB,  J ' 


fro 

u'=£l 


T4 


c, 

T4  ^ + 


7 


(3-41) 


(3-42) 


(3-43) 


and  C2  = 0,  by  the  boundary  conditions  (3-34a). 

Set  U,=f,r+  £f  , 


where 


(3-44) 


(3-45) 


/%*//  -l? 


*■/  *'■ 


Y^r/JS,  ]<* 


for  p = 1,  2,  3 


(3-46) 


With  the  help  of  the  boundary  conditions  (3- 34b),  we  can 
determine  the  unknowns , , c^  and  q^ . Consequently  Wj 

and  U}  are  obtained. 

The  same  procedure  is  applied  to  solve  for  W2  and  U2. 
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Let 


A> 


^ ^ l K 0t  7j ,/6Sa^f°'*s7^fS-'\ 

*7*0-  ■*  4* -%)ku,}^-,  i 

+[<*-*>«(% -*^+**fjg+ 

/ -d?- 


< ( <=*  f. 

A? 

=yf  ^^4 - <%. 

^S/r/ ?sf>  0,-(Sr^-,)0j ]^s/2^  - JS/cx 

(?-»[■*(**  + 00,  -seJ]°<!s°^/-*}7*'‘ 


From  Equation  (3-35a)  we  solve  for  U2  as 

&,  y f 7j  Y ^ a7  ^ ^ ? 

~/~  j5 

? • 


(3 


(3-47) 


3-48) 


-49) 
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Let  W2  be  also  expanded  into  power-series 


w*  r** 


(3-50) 


Py  substituting  of  Equations  (3-49)  and  (3-50)  into  Equa- 
tion (3-35b)  and  equating  the  coefficients  of  the  same 
order  terms,  we  get  the  following  formulae 


* ' * # c,  =(■?%  t 


(3—51 a ) 




2 /r>  B/77  (2/>7  + 2)  y.  ^ 

+0[2?O*-?+04  - 

0(f-0[2(2M-  ,3)0, 

2(2/V-2fi- 


sf-2 


4X0,  sr?(/?9-0 


+ Bj]  o'  o' 


y*  x(prj-p)  • 


for  /r?  = 2, 


-2 


(3-51b) 


Substituting  Equation  (3-50)  into  Equation  (3-49)  and 
applying  the  boundary  condition  (3- 36a),  we  have  c3  = 0 and 


<£>o 


(3-52) 
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where 


y- 


0, 


( 2^-^/rr  +-/) 


(3-53a) 


(3-53b) 


Py  applying  the  boundary  conditions  (3-36b),  these  three 


unknowns,  and  q2  can  be  determined. 

Let  us  for  a rough  approximation  assume 

l/l/,  = / 7-  ^ ^ y * j (3-54a) 

Uf  +/%73+/3r  7* , (3-54b) 

^ Y <=*£  oj  ^y  ~?/0 , (3-55a) 

^ ?7+/%79Y,,?"  . (3‘55b) 


' the 

formulae 

(3- 34b), 

(3-42) 

, (3-41),  (3-45), 

(3-46) 

d (3« 

-54)  we  obtain  the 

f irst 

set  of  algebraic 

equations 

Jl 

^2 

^ - 

_ 4 

4*,  ; 

( 3-56a) 

y - 

“/  , 

(3-56b) 

y = 

( 3-56c ) 

A 71 A ’“A 

= o 

s 

(3-56d) 

ff 

/h  - (*  = O „ 


(3-56e) 
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/j2  ) 

(s  A 63,)  ^ 

-?  — ^8, 


= O 

y 


( 3-56f ) 
(3-56g) 


Py  formulae  (3- 36b),  (3-51),  (3-53)  and  (3=55)  we  obtain 
the  second  set  of  algebraic  equations 

fj  — 3 — 

2 , (3-S7a( 


= fx*(*6, -$,)+£- (*0,t 

f£j)]oSj  - + , (3-57b) 


(3-57c) 

^ s oj;  = o , 

(3-57d) 

F + F*  i-fi. 

F /%?  + ~ c>  , 

( 3-57e ) 

-L-  f -f  0*  ] .<? 

- ^3,  ( 2/c  r ^3  y 

( 3-57f ) 

- / f . 8*  | 

S3,  ( A"  ^ , 

(3-57g) 

fir  +3  3,)  **<6 

— _Z.  r / /SB,  &s 

, -S/SB,  Bs  , *1 

(3-57h) 
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(3-571) 


3 

& 


\B?~  ~ ~ jrl  + 

* 6 '(*#  - 4(*0,  - 33, )<3?  <*]  + 

3'~  (*4  + 0j)&2°<4 


3 


&B,  K 


V ' 


F,  -(  ■ 


J-  ^ 77*J°(r 

$>  /ffj  3°°  <=3yo  = ~ /<i-3  ( 33,  ~ 3j  }cx^ ^ 

- 2-(,on  - JB,)c*J]  +JL- (64 , 3j)o£  t 


?j  , - 


')<*,  = °. 


0-57J) 

(3-57k) 


(3-571) 


A -(£-+£k)<*,.=o 


(3-57m) 


A computer  program  in  Fortran  II  (Appendix  C)  is  used  to 
solve  Equation  (3-56)  and  (3-57),  the  results  are  plotted 
as  dimensionless  loading  versus  dimensionless  central 
deflection  and  compared  with  known  solutions  (23,  24), 
as  shown  in  Figure  6. 


CHAPTER  IV 


CONCLUSIONS 

Poth  the  derivation  and  a numerical  application 
for  a large-deflection  theory  of  orthotropic  sandwich 
shallow  shells  have  been  carried  out  in  this  dissertation. 
The  numerical  solution  to  a typical  engineering  problem 
is  not  only  to  demonstrate  the  applicability  of  the 
theory  proposed  here  but  also  to  show  that  the  numerical 
procedure  adopted  is  useful  in  solving  such  sandwich 
shell  problems. 

In  the  theoretical  derivation  the  procedure  is 
general  enough  to  cover  all  the  cases  of  orthotropic 
sandwich  shallow  shell  structures.  All  the  quantities 
involved  are  expressed  in  general  coordinates  so  that 
one  can  easily  apply  them  to  any  specific  geometrical 
configuration  by  the  proper  coordinate  transformation. 

The  derivation  is  based  on  the  following  considerations. 

(1)  The  structure  is  shallow;  i.e.,  the  rise 
of  its  middle  surface,  z,  and  its  first  as  well  as 
second  derivatives  are  so  small  that  their  products  are 
negligible  in  comparison  to  other  dimensions. 
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2)  The  facings  of  the  structure  are  thin  enough 
that  Kirchhof f-Love  hypothesis  is  valid. 

7)  he  transverse  shear  deformation  affect  is 
included  in  the  core. 

4)  The  planes  normal  to  the  middle  surface  re- 
frain plane  but  no  longer  normal  to  the  middle  surface 
after  deformation. 

5)  No  transverse  compression  through  the  thickness 
is  considered,  i.e.,  the  deflection  w is  independent  of 
the  transverse  coordinate. 

6)  A generalized  plane-stress  theory  is  applied. 

7)  The  temperature  distribution  is  a known  func- 
tion of  the  space  coordinates. 

Most  of  these  considerations  are  commonly  used  in  ordinary 
shell  theory  except  for  items  (3)  and  (4)  which  distinguish 
sandwich  shell  theory  from  classical  shell  theory. 

For  the  numerical  solution,  the  temperature  terms 
are  dropped  for  simplicity.  However,  this  does  not  cause 
any  difference  in  the  numerical  method  to  be  used  since 
the  temperature  is  a known  function  which  effects  only  a 
finite  number  of  terms  in  the  series  solution  and  the 
recurrence  formulae  still  hold. 

No  effort  has  been  made  to  improve  the  accuracy 
of  the  numerical  result,  since  this  can  be  done  easily 
by  increasing  the  number  of  terms  in  the  expansion.  Even 
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for  a rough  approximation,  the  result  is  within  10  per 
cent  of  known  solutions.  The  problem  is  completely 
solved  once  the  deflection  w and  the  displacement  u and 
the  angular  rotation  0 are  determined. 
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TABLE  1 

COEFFICIENTS  FOR  SANDWICH  SHALLOW  SPHERICAL  SHELL 

k = 2.0 


Oio 

1.0 

/*/ 

0.53 

y, 

5.23 

^*2 

-2.70 

P* 

-1.82 

n 

-10.3 

2.55 

Ps 

1.99 

?s 

5.19 

^6 

-0.85 

P? 

-0.71 

y7 

-0.03 

C*o 

0.0 

0.24 

y, 

-2.15 

1.16 

p3 

-0.44 

% 

10.39 

-2.81 

Ps 

0.49 

ys 

-11.37 

2.05 

Pr 

-1.05 

^7 

0.50 

-0.17 

P* 

1.25 

y9 

2.68 

-0.24 

A 

-0.50 

% 

-0.05 

90 


TABLE  2 

COEFFICIENTS  FOR  SANDWICH  SHALLOW  SPHERICAL  SHELL 

k = 1.5 


cto 

1.0 

P> 

0.53 

y, 

5.31 

-2.74 

-1.84 

y, 

-10.56 

cJ4 

2.61 

2.04 

>s 

5.27 

-0.87 

A 

-0.72 

^7 

-0.02 

0.0 

0.29 

y, 

-1.63 

0.88 

ft 

-1.09 

y, 

8.08 

-2.19 

2.02 

7* 

-9.18 

1.66 

-2.47 

y7 

"0.78 

-0.18 

F 

1.81 

y? 

1.98 

^,0 

-0.17 

F" 

-0.56 

r„ 

-0.04 

91 


TABLE  3 


COEFFICIENTS  FOR  SANDWICH  SHALLOW  SPHERICAL  SHELL 

k = 1.0 


Oio 

1.0 

0o  53 

y, 

5.37 

<*2 

=2.77 

ft 

-1.86 

-10.67 

2.65 

ft 

2.07 

$ 

5.32 

“0.38 

A 

-0.74 

>T 

-o.oi 

o<D 

0.0 

0.40 

y 

-1.04 

ft 

0.57 

ft 

-2.20 

% 

5.42 

-1.47 

ft 

4.59 

£ 

-6.44 

ft 

1.17 

ft 

4.59 

% 

0.78 

ft 

-0. 15 

ft 

2.91 

% 

1.30 

<^JO 

-0.12 

P» 

-0.72 

% 

-0.03 
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TABLE  4 

COEFFICIENTS  FOR  SANDWICH  SHALLOW  SPHERICAL  SHELL 

k = 0.5 


O<0 

1.0 

ft 

0.52 

5.40 

-2.79 

ft 

-1.87 

>3 

-10.75 

*4. 

2.68 

ft 

2.09 

>r 

5.35 

ft 

-0.89 

A 

-0.74 

^7 

-0.002 

C*0 

0.0 

A 

0.76 

y 

—0.41 

o!^ 

0.22 

A 

-5.08 

% 

2.51 

ft 

-0.69 

Pr 

11116 

£ 

-3.32 

0.60 

ft 

-11.65 

% 

0.57 

C/q 

-0.08 

A 

6.14 

0.66 

—0.06 

A 

-1.32 

% 

-o.oi 
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M SUREACE 


FIG.  I PROJECTION  OF  M-SURFACE  ON  TT-PLANE 
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P1 


FIG.  2 DISPLACEMENT  VECTORS 
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FIG-  3 


SANDWICH  STRUCTURE 
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W 
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FIG.  4 SANDWICH  SHALLOW  SPHERICAL  SHELL 
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FIG.  5 

VARIATION  OF  CENTRAL  DEFLECTION 
WITH  LOADING,  SANDWICH  SHALLOW 
SPHERICAL  SHELL 


1.0 


1.5 


2.0 


NON-DIMENSIONAL  CENTRAL  DEFLECTION  vvo  = (w/h)rr0 


35 
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APPENDICES 


APPENDIX  A 


APPENDIX  A 


PROOF  OF  THE  CONVERGENCE  OF  THE  POWER  SERIES 

From  the  recurrence  formula  (3=18),  solving 

for 


/ 


L 


4 


ex'. 


-2 


(A-l) 


where  , b2  and  b3  are  defined  by  Equation  (3-20n). 
Since  % is  in  the  interval  of  [o , l]  , let  us  only  con 
sider  the  ratio  test  of  c*8s. 

Replacing  p by  p + 1 and  by  substitution  of 
Equation  (A-l),  we  have  the  coefficient  for  the  next 
term 


Of 


(if  ' A,) 


of. 


-+ 


__  h 

27r?+3f(?^y7  ^-2  . 

r ' (A-2) 
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The  ratio  is 


= /—■ 7 - /■  (ft  + ^)a  f 
1 2C?*jfJ  / ' 4 * y a*?>-  - 


; tz 


i±-£> 


+ JLEJL7-:  ce  - 


o' 


4 _ 4 


C*^,  + 


Of„ 


ytz  • jy„/y  %'*) , (A.3) 


There  are  two  possibilities: 


Case  1 


ex' 


-s 


s' 


o(. 


Dividing  both  the  numerator  and  denominator  by 


2p— 2*  we  have  the  ratio 


■{«*¥) 


Jbf  2 


CX 


+ 


2^-2 


*$lk 

i 

Nl 

*>  / / 

**(?+*/?*  b}/ 

f /*>?■? 

, ** 

4 

°<2f-2  &(?«/?*. 

(A-4) 


Let  us  take  the  limit  of  Equation  (A-4)  as  p becomes  very 
large.  Suppose  ^Q^,^/^<a^^does  not  approach  to  zero 
as  p— ■ Otherwise  the  series  is  convergent  by  this 

as  p -rc>o 


condition.  Since 
C^s/b  ts 


implies 


ex. 


<^2j6t2 

_ 

cXjA+S 

*** 

^Sfb 

i 

*°  < / ✓ 


as  p -r-  0o  . 
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JL, 

fi 


o( 


^ t 6 


O 


■?fi*2- 


fi~>&° 


2(?tsf 


' /im 


't '££**-  ' 


(—  - ^ ) rv 

( h 2 / CV- 


Z_  _ 


h 


&(?+0*  ^fi-Z  ~ 32(fi+tf?% 


( 


& 2 /it  2 


4 


cy 


h_ 


°V-J  o/^.2 


=2„ 
fi  -*t>o 


/ 

±+  *2 

2 Cfit3/ 

4 ; + 

(A-5) 


Since  b^ , b2  and  b3  are  fixed  constants  which  are  indepen- 
dent of  p,  at  a large  value  of  p we  can  have 


o( 


^ fi  t 6 


<y 


Jfit* 


o 


/ 


( A— 6 ) 


Hence,  in  this  case,  the  series  is  convergent  by  the  ratio 
test . 


Case  2 


o( 


2fit  2 


o( 


Let  us  divide  the  numerator  and  denominator  by 


2p  + 2 


J 


cy 


Jfi  + 4 *(fit3J* 

- >y 


fi*!,  CXjpn  j/ / / + g(fi„f 


**  <=*££ 


c2 


32  ( & 

y-2 


■} 


( A-7 ) 
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Taking  the  limit  of  its  absolute  value  as  p becomes  very 
large,  we  have 


Cx* 


/■  6 
+4 

•f  ~ 


-¥) 


?(  p+/j‘ 

k* 


f'+ 


ft  -*  ^ 


/ 

(a  *&+ 

f 

kj 

/ 

' ^ 4 

4 * 

(A-8) 


Py  the  same  argument,  we  can  see  that  the  series  is  con- 
vergent also  for  Case  2. 


APPENDIX  B 


COMPUTER  PROGRAM  2 


SANDWICH  SHALLOW  SPHERICAL  SHELL 

DIMENSION  A ( 25 , bO ) , B(25,50),  C(25,25),  F(25),  G(15), 

H { 25 ) ,X( 10) 

1 READ  INPUT  TAPE  5,2,M,(X(I),  1=1, M) 

2 EORMA  T ( 115,  9F8.3) 

WRITE  OUTPUT  TAPE  6,3,M,  (XII),  1=1, M) 

3 FORMAT  ( /// 40X , 5H  DATA//40X,  1 4/ ( 36X , E 1 A . 4 ) ) 
BA=(2.*X(3)-4.»X(i) ) *X ( 8 ) *X { 8 ) 

BB=2 . * ( X ( 1 ) +X  { 2 ) -X  ( 3 ) ) *X ( 8 ) *X { 8 ) 

8C=X ( 4 ) +X ( 6 ) 

BD=X ( 5 ) +X ( 6 ) 

B6=X ( 4 ) /X ( 7 ) -X ( 6 ) *BC/X  ( 7) 

BF  = BA*X(3)*X(4)/(XU)*X(  7)  ) 

BQ=(0.5*XI 1 ) +.25*X( 3 ) ) / X ( 8 ) 

YA=2.*X{ 1 ) -X ( 3) 

YB=4.*X( 1 ) -X ( 3) 

YC=YB+2. *X{ 1 ) 

YD=YC+2.*X( 1 ) 

YE= 10 . *X ( 1 )-3.*X( 3) 

YF=BA*X(4)/(BC*X( 7) ) 

YG=X ( 6 ) - 1 . 

YH=BA*X( 4) / ( X ( 1 ) *X (7 ) *X ( 8 ) ) 

C SET  UP  FIRST  MATRIX 

N=  1 3 
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00  10  1=1, N 
DO  10  J = 1 , N 
10  A { I , J )=0 
A ( 1 , l ) =0 . 

A ( l,2)=-2.*BF 
A ( 1,3) =-32 • * C l.*BC) 

A { 1 »4 ) =1 152 • *BE 
A ( i, 9)=,5*BA 
A ( 2 » 2 ) = 1 . 

A( 2 ,3)=1 . 

A ( 2 , 4 ) = 1 . 

A ( 3,2)=.25«BF 
A { 3 , 3 ) =4 . *BC 
A(3,4)=-L44.*BE 
A{  3, 11 )=BC 

A ( 4 » 2 ) =*2  • * YG*  ( 1 • *X  ( 9 ) ) 

A( 4 » 3)=4.*YG*( 3 • +X ( 9) ) 
A(4,4)=6.*YG*(5.+X(9) ) 

A( 4 , 10 )=X ( 6 ) 

A( 4, l L ) =3.*X( 6) 

A(4, 12)=5.*X(6) 
A(4,13)=7.*X(6) 

A(  5 , 5 ) 3 L • 

A ( 5 , 6 ) = 1 . 

A ( 5 , 7 ) =1 . 

A ( 5 , 8 ) = l . 

A(6,2)=-{0.5+0.25*X(3)/X( l ) J 
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A ( 6 , 6 ) = 1 . 

AC7,3)=-{2./3.  +X{ 3 ) /( 6. *X  ( 1 ) )) 
A ( 7 » 7 ) = l . 

A ( 8 » A ) =- ( .75+.125*X(3)/X< 1)  ) 

A ( 8 » 8 ) =1 . 

A(9,5)=l. 

A( 9,9) =-0.5 
A ( 10,  1 )=.5*YF/BA 
A { 1 0, 2 ) =2 . 

A ( 10, 3 )=-32.*BL'/RC 
A( 10,9)=.5*YF 
A ( 10,10)=1. 

A( 11,10 )=1. 

A ( ll,li)  = l. 

A ( ll,12)=l. 

A ( ll,13)=l. 

A ( 12»3)=BF/(6.*BC) 

A ( 12,4) =6 . 

A ( 12,12)=!. 

A( 13,4)=. 125*BF/BC 
A( 13,13)=!. 

CALL  INVtRT  (N,A,C) 

DO  19  1=1, N 
L5  F ( I ) =0 

F ( 1 )=-.5*BA*X{ 3)/X( 1 ) 

F { 2 )=- 1 . 


F(9)=.5*X{3)/X( 1 ) 


H 10)=-.5*BF/BC 
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0G  20  I = l , N 
SUMA=0. 

DO  19  J = 1,N 

19  SUMA=SUMA+C ( I y J 1 *F ( J 1 

20  G( I )=SUMA 

WRITE  OUTPUT  TAPE  6V4«N«(G(I)»  1=1, N) 

4 FORMAT  ( / / /45X , 5HF IRST1X » 3HSET1X » 12HC0EFFICIENTS//10X# 

I5/( 20X.4E20 

I. 91) 

C SET  UP  SECOND  MATRIX 

N=  1 9 

00  35  1=1 ,N 
DO  35  J=1 ,N 
35  A ( I » J ) =0 
A { l, 11=0.5 
A I 1 # 2 1 =-2 . *BF 
A ( 1 , 3 1 =-32 • » ( 1 • +BC 1 
A ( ly 41=1152. *BE 
A ( 1 ,9)=.5*BA 
A(2,2)=l. 

A ( 2 , 3 1 = 1 . 

A { 2 » 4 1 = 1 . 

A(2yl4)=i. 

A(2, 15)=L. 

A( 3,2)=.25*BF 
A ( 3,3)=4.*BC 


A { 3 y 4 1 =-144 . »BE 


AC  3» 11 )=BC 
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AC4,2)=2.*YG*(  l • ♦ X { 9 ) ) 
A(4,3)=4.«YG*(3.+XC9) ) 
AC4f4)=6.*YG*{5.+X(9)  ) 

A ( 4 * 1 0 ) =X { 6 ) 

AC4»11)=3.*X{6) 

A(4» 12 ) ~5. *X ( 6 ) 

A{  4 » 13 ) =7 • *X  ( 6 ) 

A(4»14)=a.*YG#( 7 . +X ( 9 ) ) 
A(4»L5)=10. *YG* ( 9. +X ( 9 I ) 
A(4t18)=9.*X(6) 

A(4tl9)=ll.*X(6) 

A ( 5 , 5 ) = L . 

A ( 5 , 6 ) = 1 . 

A { 5 *7 ) = 1 . 

A ( 5 , 8 ) = 1 . 

AC  5» 16 ) -1 . 

A (5, 17) -1. 

A(6,2)=-(0.5+0.25*X{  3 ) / X C 1 ) ) 

A ( 6 f 6 ) = 1 . 

AC  7,3)*- 12./ 3.  +X(3)/(6.*X(  1)  ) ) 
A ( 7 , 7 ) = i . 

AC  8, 4 ) =— { .75  + . 125* X( 3) /X(  1 ) ) 

A ( 8 » 8 ) = 1 . 

AC9,5)=1. 

AC9f9)=-0.5 


AC  10*1 ) = .5*YF/BA 


A ( 10, 2 ) =2 
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A{  10,3)=-32.*BE/BC 
Al 10»9)=.5*YF 
A( 10, L 0 ) = 1 . 

A ( 11,  I0)  = l. 

A ( 1 1 , 1 1 ) = 1 . 

A ( 11, 12 1*1. 

A(ll,13)=l. 

A ( 11, 18)=1. 

A { l 1,  191*1. 

A( 12,3)=BF/(6.*8C) 

A { 12,4) =6 . 

A ( 12, 12 ) =1. 

A ( 12, 14)=-384.*3E/BC 
A ( 13,4)=.125*BF/BC 
A ( 13, 131=1. 

A { 13,14)=8. 

A{ 13, 15)=-800.»BE/BC 
A ( 14,2)=.25*BA*X{3)/X{  1) 

A ( 14,3) = -4 . * BF 
A(  14,4) *■— 1 44 . * { l.  + BC) 

A ( 14, 14 )=9216.*BE 
A ( 15,3)-BA*X(3)/(6.*X{ 1 ) ) 

A { 13,4) =-6 . *BF 
A ( 15, 14) =-384. *< 1.+3C) 

A ( 15,  15) =38400. »BE 

A( 1 6 , 1 4 ) =-  ( .8*.1«X(3)/X{  1)  ) 


A ( 16, 16 ) = 1 . 
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A( 17, l5)=-(5./6.+X(3)/( I2.*X(1) ) ) 

A( 17, 17)=l. 

A ( 18, 14 ) 3 • 1 *GF/BC 
A ( 18, 1 8 ) = 1 0 . 

A ( 18, 1 8 ) =1 . 

A ( 19*15)=BF/{ 12.*8C) 

A(19,19)=1. 

CALL  INVERT  (N,A,C) 

DU  AO  I = 1 , N 
AO  F ( 1 ) =0 

F ( l ) =4. *X{8)*X(4)«YA*(8.*G(b)*G(3)+A.*G(2)*G{6)-3.*G(2 

) *0 ( 2 ) )/ 

lX(7)-2.*X(8)*YA*G(5)*G(2)-2. *YH« ( X { 1 ) + . 5*X (4 ) ) *G( 2 ) »G( 

2) 

F(3)=*-X(  A)»X(8)*YA»(A.*G(5)*G(3)  +2  . *G  ( 2 ) »G  ( 6)-l . b*G  ( 2 ) 

*G  1 2 ) )/X( 7) 

1 + . 5*YH«  ( . 8*  X ( 1)+«25*X(3I  ) *G(  2 ) *G  ( 2 ) 
F(6)=-.5*6Q*G(2)«GC2)/X(1) 

F(7)=-(2.*X(l)+.5*X(3))*G(2)*G(3)/(3.*X(l)*X(8)) 

F ( 8 ) =- ( ( .75*X(1)  + .125*X(3))*G(2)*G(4)  + (.5*XC1)+X (31/12 

. ) *G ( 3 ) * 

1G(  3 ) ) / ( X { 1 ) * X ( 8 ) ) 

F( lO)=-2. »YF*YA*X( 8)*G(5)*G(2)/CA 

F( 12)=-YF*X(8)*(6.*YA»G(5)*G(4)+(4.*(YD-2.*X( 3) )*G(3)* 

G ( 6 ) +6 . * YA 

l*G( 2)*G( 7)-4.*( YD- 2. *X ( 31 )*G(2)»G13) 1/3. I/BA+YH* (2.*X( 

1 ) + • 5 *X ( 3 ) ) 

2*G(3)*G(2)/(3.*BC) 

F(  13)=-.25*YF*  ( ( 3.*X( 1 ) — X { 3 ) ) *(24.*G(4) *G(6) +16.*G( 3)* 

G( 7) ) +8 . * 
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1YA*G<2)*G(8 )-YE*(4.*Gl3) *G(3 )+6.*Gt 2) *G(4 ) ) J/BA  + YHM3. 

»Xtl ) + .5* 

2X ( 3 ) ) * ( .25*G12)*G( 4) +G t 3 ) *G ( 3 ) /6 . ) /8C 

F( 14  ) =X  t 4 ) *X  ( 8 ) «•  ( YA* ( 144.*G(5)*G14) +48.*G(2) *G{7) ) *32. 

*(  ( 8 . *X ( 1) 

1-3. *Xt 3) )*G(3)*G(6)-(22.*X( 1) -9. *X(3))*Gt2)*G(3) ) ) / X ( 7 

) -X { 8 ) * Y A* ( 

24. *G( 5 ) *Gt  3 ) +2 . *Gt  2 ) *G l 6 ) -1 • 5*G( 2 ) *G l 2 ) J -4. *YH*  1 4. *X ( 1 

) ♦X ( 3) )* 

3G(2)*G(3)*.5«BA*3Q«G(2)*G(2)/X( 1 ) 

F{ 15) =24. *Xt8)*X (4 )*( ( 3 . *X ( i ) -X ( 3 J ) * ( 12 . *G ( 4 ) *G( 6 ) +8. * 

G( 3) *G( 7 ) ) + 

14. * YA*G l 2 ) *G ( 8 ) -6. *YE»G  C 3 ) *G  1 3)  — 12 • * ( YO-2 . *X ( 3 ) ]*Gt2)* 

G { 4 ) ) /X { 7 ) — 

26 . * YA* X 1 8 ) *G 1 5 } *G ( 4 ) - X t 8 ) * C 4 . * { YD- 2 • *X ( 3 ) )*(Gt3)*G(6)- 

G ( 2 ) *G( 3 ) > + 

36.*YA*G(  2 ) *G(  7)  ) /3.-YH*  ( 6.*X(  l)  + X(  3)  )*  (4.  *G(  3)  *G(3  J +6. 

*G ( 2 ) *G ( 4 ) ) 

4+8A*(4.*X( l)+X(3))*G(2)«Gt3)/(6.*X(8)»Xtl)) 

Ft 16)=-3.*(8.*X( 1)+X<3) )*Gt 3)*G(4)/(20.*X(8)*X(  1) } 

Ft  1 7 ) -- 3 . * ( 10.*X(1)+X(3) )»Gt4)*G(4)/l40.*X(8)*X(l) ) 

Ft  1 8 ) =— X ( 4) *X ( 8 ) * ( 6. * t YD+YE-X 1 3 ) ) »G 1 41 *Gt  7 ) +4.*( YE+YC- 

X ( 3 ) ) * G t 3 ) 

l*Gt  8J-36. *YB*Gt  3 ) *Gl 4 ) ) / { 5 . *X ( 7) «BC ) + • 3« YH* 1 4. *X ( 1 ) + .5 

* X ( 3 ) ) * G ( 3 ) 

2»G14)/BC 

Ft  19 )=-X ( 4 ) *X ( 8 ) * ( 12.*YB*Gt4) *G t 8 )-3 . * l YC+YD-X ( 3 ) )*Gt4 

) *G ( 4 ) ) 

l/t2.*Xl7)*BC)*.3*YFI*(2.5*Xt  l ) *.25*X  t 3)  ) *G  ( 4)  *G  (4 ) /BC 
DO  50  K=1 ,N 
SUM=D 

DO  41  L = l , N 

41  SUM  = SUM  *-C  t K , L ) *F  { L ) 


50  H{  K ) = SUM 
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WRITE  OUTPUT  TAPE  6 , 5 , N , (H (I) , I =1 , N) 

5 FORMAT! ///44X, 6HSEC0ND1X , 3HSET1X, 12HC0EFFICIENTS//10X, 

1 5/ ( 20X , 4E20 

1 « 8 ) //  ) 

wo=o 

DU  60  I = L , 10 
WQ=W0*0.2 
FAC  T=l./128. 

Q=F AC T * ( G( 1 ) +H(  L ) « WO  J * WO 
60  WRITE  OUTPUT  TAPE  6,6,W0,Q 

6 FORMAT ( 30X,2E25. 8 ) 

GO  TO  1 

END 

SUBROUTINE  INVERT  (N,A,C) 

DIMENSION  AI25.50),  B(25,50),  CI25.25) 

K = N 

KF  = N + 1 

KK  = N + N 

DO  30  I = 1,  N 

K * K «-  I 


DO 

31 

31 

A { I 

f J ) 

30 

All 

, K) 

DU 

32  I = 

It 

32 

8 ( I 

, l)  = 

A ( I 

R = L.O/At 1,  1 ) 


DO  33  I 


2,  KK 
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33  B(  1,  I)  = R « A(  1,  I ) 

DO  34  I = 2,  N 

DO  35  J = I » N 
B(  J,  I ) = A(  J,  I ) 

LAST  = I-  l 

DO  36  Jl  = 1,  LAST 

36  B ( J , 1 ) = B ( J , I ) - B(JyJl) ) 

35  CONTINUE 

KUL  =1+1 

DO  37  J = KOL,  KK 

B(  I,  J)  * A(I,  J) 

DO  38  Jl  = 1,  LAST 

38  8(  1 , J)  = B l I » J ) - BIJl,  J ) * ft  ( l , J l ) 

B( I, J)  =B( I , J)/Bl  1,1  ) 

37  CONTINUE 

34  CONTINUE 

K1  = N+  L 
NT  = N - I 

DO  39  1 = 1,N 

* 

NK  = N *■  I 
C ( N , I ) = B ( N,NK ) 

DO  40  J = 1 , NT 

NJ  = N - J 

CINJ,  I ) = B ( N J , NK ) 

DO  41  J3  = 1,  J 
N JL  = K1  - J3 


41 


CINJ,  I) 


CINJ,  l) 


CINJL,  I)*B(NJ,  NJL ) 
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40  CONTINUE 

39  CONTINUE 
RETURN 
END 


DATA 
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COMPUTER  PRUGRAM  1 
SHALLOW  SPHERICAL  SHELL 

DIMENSION  A ( 25*  50  ) » B(25,50),  C(25,25),  F{25),  GU5), 

H T 25 ) ,X(lO) 

1 READ  INPUT  TAPE  5,2,M,(X(I),  I=I,M) 

2 FORMAT  (1110,  8F8.3) 

WRITE  OUTPUT  TAPE  6,3,M,  (X(I),  1=1, M) 

3 FORMAT  ( /// 40X  f 5H  DAIA//40X,  1 4/ { 36X , E 1 4.4 ) ) 

BA=(2.*X( 3)-4.*X{ I) ) *X ( 8 ) *X { 8 ) 

BB  = 2 . * ( X { 1 ) +X ( 2 ) -X ( 3) ) » X ( 8 ) »X ( 8 ) 

BC=X  ( 4 ) +X ( 6 ) 

BD=  X ( 5 ) +X ( 6 ) 

BE=X(4)/X(7)-X(6)*BC/X( 7) 

GF=8A*X ( 3 ) *X ( 4 ) / { X ( I ) * X { 7 ) ) 

BQ= ( 0. 5 * X ( 1 ) + • 25»X ( 3) )/X{8) 

YA=2 . *X ( 1 )— X ( 3 ) 

YB  = 4 . * X { 1 ) -X  ( 3 ) 

YC=YB+2.*X( I ) 

YD=  YC+2 . *X ( l ) 

YE=10.*X(1)-3.*X(3) 

C SET  UP  FIRST  MATRIX 

N=7 

DO  10  I = 1 » N 
DO  10  J = 1 , N 


10  A { I » J ) =0 


A ( t,  1)  = .5 
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A ( 1 » 3 ) =-32. 

A ( 1, 7)=.5*BA 
A ( 2 » 2 ) = 1 . 

A(2,3)=l. 

A ( 3 » 2 ) = 1 . 

A ( 3 » 3 ) = 2 • 

A(4,4)=i. 

A ( 4 f 5 ) = l • 

A ( 4 , 6 ) = 1 . 

A(5,5)=l. 

A(5,2)=-( .5+.25*X(3)/X(i) ) 
A(6,6)=l. 

Al6,3)»-l2./3.tX(3)/{6.*X{l))) 
A( 7, 7)=-. 5 
A ( 7 * 4 ) = l . 

CALL  INVERT  IN,A,C) 

DO  15  1=1, N 
15  F ( I ) =0 

Ft 1)--.5*BA*X(3)/X( 1 ) 

F { 2 ) =-  l . 

Ft  7 ) = .5*X( 3 ) /X( i ) 

DU  20  1=1, N 
SUMA=0. 

DO  19  J = 1 , N 

19  SUMA=SUMA+C( l,J)*F(J) 


20  G ( 1 ) =SUMA 
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WRITE  OUTPUT  TAPE  6»4,N,(G(l),  1=1, N) 

4 FORMAT  ( ///45X » 5HF IRST IX, 3HSET1X, 12HCOEFF1C IENTS// IOX, 

I 5/ ( 2GX ,4E20 

1.8)  ) 

G ( 9 ) =G  { 7) 

G(8)=Q 
G( 7 )=G( 6 ) 

G(6)=G< 5) 

G( 5 )=&( 4) 

G ( 4 ) =0 
G( 10)=0 
G( 1 1 ) =0 

C SET  UP  SECOND  MATRIX 

N=  1 3 

DU  35  1=1, N 
DO  35  J = 1 » N 
35  A ( I » J ) =0 
A C i , l)=0.5 
A ( l,?)=-2.*BF 
A ( 1 , 3 ) =-32 . * ( 1 • +BC ) 

A ( 1,4) =1152. *BE 
A ( 1 » 9 ) = . 5 *BA 
A { 5 , 5 ) = 1 . 

A { 5 , 6 ) = 1 . 

A ( 5 , 7 ) = 1 . 

A ( 5 , 8 ) = 1 . 

A ( 6 , 6 ) = 1 . 


A(6,2)=-{0.5+0.25*X( 3)/X{ 1 ) ) 
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A ( 7 « 7 ) * l . 

A(7»3)=-(2./3.  *X< 3)/(6.»Xtl) ) ) 
A(8»8)=l. 

A(8,4)=-t.75+.125*Xt3)/Xtl) ) 

A ( 9 » 5 ) = 1 . 

A { 9 » 9 ) =-0. 5 

A{2,2)=.25*BA*X13>/X« 1 ) 

A{  2, A )=- 144. 

At  3* 2 } = 1 . 

A l 3 » 3 ) = 1 . 

A ( 3 , 4 ) = 1 . 

A(3,I0)-1. 

A ( 3 » 1 2 ) - 1 . 

At4,2)=l. 

A ( 4 , 3 ) = 2 • 

A(4,4)=3. 

A(4,10)=4. 

A ( 4 , L 2 ) = f> . 

A(5»il)=l. 

A(5, 13)*1. 

At 10,3)=8A*X(3)/(6.*X{  i)  ) 

At  10, 10 ) =-384 • 

At  11, 10 i =- t . 8+.i«X(3)/X( 1) ) 

At  11,111  = 1. 

At  12,12) =-800. 

At 1 2 , 4 ) = . 1 2 5 * B A * X ( 3 ) / X ( 1) 


At 13,l2)=-{5./6.+X(3)/(12.»Xti)) ) 
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A ( 13*13 )=i. 

CALL  INVERT  (N,A,C) 

DU  40  1=1, N 
40  H I ) =0 

F( l )=4.*X( 8) *X{ 4) *( B.*YA*G( 5) *G( 3 ) +4 . » YA*G ( 2 ) *G( 6 ) -3. * 

YA  *G ( 2 1 

l*G(2) )/X(  7)-2.*X(8)*YA*G(5)*G(2)+2.*BA*X(4)* (X(1)+.5*X 

{ 4 ) ) *G ( 2 ) 

2*GI 2) / (X { 1 ) *X ( 7 ) *X ( 8 ) ) 

F(2)=( 1.5*XI8)*YAfBA«(2.*X( 1)+X(3) ) / ( X ( l) *8 . »X { 8 ) ) ) *G( 

2 ) *G( 2 )- 

1X(8)*(4.*YA*G(5)*G(3) *2 . * YA *G ( 2 ) *G ( 6 ) ) 
F(6)=BQ*G(2)*G(2)/{2.*X( 1) ) 

F ( 6 )=-F ( 6 ) 

F ( 7 ) = ( ( 8Q+X  { 1 )/X(8)  )*G(2)*G(3H-BQ*G(2)*G(3)  )/(3.*X(  1 ) ) 
F ( 7 ) =-F ( 7 ) 

F ( 8 ) = ( •25*(&Q+2.*X(1)/X(6) ) *G (2 ) *G (4 ) + ( BQ+X ( 1 ) /X ( 8 ) ) *G 

(3)»G(3)/3. 

l + . 25*3Q*G( 2 ) *G ( 4 ) ) /X ( 1 ) 

F ( 8 ) =-F ( 8 ) 


F ( 10)=-6.*X{8)*YA»G{5)*G(4) ♦X(8)»(4.*(8.*X(1)-3.*X(3)) 

*G ( 3 ) *G(6)+ 

16. «YA*G( 2 ) *G( 71-4. *( 8.*X( 1 )-3.*X(3) )*G(2) *G(3) 1/3.+BA* 

{ 4 . *X ( 1 ) 


2*X( 3) ) *G( 3 ) *Gl 2 ) /« 6.«X( l)*X(8l) 

F ( 12 )=-X  C 8 1 * ( ( 12.*X( l ) -4. *X ( 3) )*G(3)*G( 7>-1.5»( 10.*X(1 

) -3. »X I 3 ) ) 

l*G( 3)*G( 3) ) + BA* ( 6. *X ( 1 ) +X ( 3 ) ) *G( 3 ) *G C 3 ) / ( 12. *X ( l ) *X ( 8 ) 

) 


DO  50  K = 1,N 


SUM-0 


124 

DU  41  L= 1 » N 
41  SUM=SUM+C ( K, L ) *F ( L ) 

50  H { K ) =SUH 

WRITE  OUTPUT  TAPE  6, 5, N, ( H (I) , 1=1 , N) 

5 FORMAT ( ///44X,6HSECUND1X,3HSET IX, 1 2HCUE FF I C 1 EN IS // IOX , 

I 5/ ( 20X,4C20 

1.8)//) 

wo=o 

DO  60  1=1,10 
W0= WO+O . 2 
FAC  T= l . /64 . 

Q=  FAC  TM  G(  L ) +U(  1)*W0)*WU 
60  WRITE  OUTPUT  TAPE  6,6,W0,Q 

6 FORMAT (30X,2E25,8) 

GO  TO  1 

END 

SUBROUTINE  INVERT  ( N » A , C ) 

DIMENSION  A ( 25 » 50 ) , 6(25,50),  C(25,25) 

K = N 

KF  = N + 1 
KK  = N N 
DO  30  I = 1,  N 
K = K + 1 


DO  31 

3 L 

All,  J) 

30 

A(  I,  K ) 

DO  32  I = l,  N 


32  8(1,  l ) 


A ( 1,1) 
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R = 

1. 

0/A { 1, 

l) 

DO 

33 

1 = 2 

, KK 

33  B(  l 

. 1 ) 

= R * 

A{ 1,  I ) 

DO 

34 

I = 2y 

N 

DO  35 

J = I • N 

B ( J , I ) = A ( J , 

LAST  = 

I-  l 

DU  36  J1  = 1,  LAST 

36  B ( J • I ) = B( J i [ ) - B(J,J1)*B( Jl,I) 

35  CONTINUE 

KUL  =1+1 

DO  37  J = KOL,  KK 

B(  I,  J)  = Ad,  J) 

DO  38  J1  = 1,  LAST 

38  0 ( 1 ( J ) = B ( I f J ) - B(J1*  J ) *B ( I » J L ) 

B ( 1 ( J ) =B ( 1 » J ) /B ( I » I ) 

37  CONTINUE 
34  CONTINUE 

Ki  = N+  l 
NT  = N - 1 
00  39  I = 1 ,N 
NK  = N + I 
C ( N » I ) = B ( Nf  NK ) 

DO  40  J * 1 » NT 

NJ  = N - J 

C(NJ,  1 ) = B (NJ,NK) 

DO  41  J3  = It  J 
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NJL  a KL  -J3 

41  C ( NJ • I)  ^ C(MJ,  I)  - C(NJl,  I ) *B ( N J » NJL) 

40  CONTINUE 

39  CONTINUE 
RETURN 
END 


L)A  1 A 
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